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ABSTRACT 


This thesis is a continuation of optimization modeling research conducted at the Naval 
Postgraduate School for the U.S. Air Force Studies and Analyses Agency. That work 
resulted in Throughput II, a multi-period model for determining the maximum on-time 
throughput of cargo and passengers that can be transported with a given fleet over a given 
network, subject to appropriate physical and policy constraints. This and other existing 
deterministic strategic airlift models assume all data are known prior to making a decision; 
often times these assumptions are unrealistic. One such assumption is aircraft reliability. 
This thesis addresses the uncertainty of aircraft reliability, which, if ignored, can result in 
models that are overly optimistic with respect to throughput capability. To address this 
issue, this thesis adds a stochastic eieiiiees to Throughput II, resulting in a two-stage 
stochastic linear program with recourse that is solved using Benders decomposition. To 
analyze the stochastic program, a simulation model of the strategic airlift system is also 
developed. This simulation model allows the user to analyze the deterministic and 
stochastic models and to compare solutions. The stochastic model, in addition to the 
features of Throughput II, accomplishes the following: (1) selection of aircraft routes by 
anticipating potential bottlenecks in the system, (2) prevents unreliable aircraft from using 
capacity limited airfields and (3) a flow of cargo from origin to destination that 1s not 


interrupted by the random events of aircraft reliability. 
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EXECUTIVE SUMMARY 


The magnitude of the airlift effort during Desert Shield/Storm was 
unprecedented. At the height of the war during Fall 1990, the Air Force averaged 17 
million ston-miles per day of cargo and troops and by 10 March 1991, strategic airlift had 
moved more than 500,000 people and 540,000 stons of cargo (Gulf War Air Power 
Survey, 1993, p. 3). If the United States was to experience a future Operation Desert 
Shield/Storm type scenario, massive amounts of equipment and large numbers of 
personnel would have to be transported over continents and oceans, with an impending 
deadline. The magnitude of such a deployment imposes great strains on current and 
future air, land and sea mobility systems. 

After the Gulf War in 1991, Patel commissioned a Mobility Requirements 
Study (MRS) of the United States Armed Forces. The MRS studied all aspects of 
mobility (domestic transportation, inter-theater lift, intra-theater lift and prepositioning) 
to determine the proper mix of sea, air, and amphibious lift, surface transportation and 
prepositioning. The goal was to provide Congress with an integrated plan for procuring 
the necessary lift for power projection in the 21st century. Two linear programming (LP) 
optimizations models were developed as part of the MRS. They are the Mobility 
Optimization Model (MOM) developed by the Joint Staff's Force Structure, Resource 
and Assessment Directorate (J8) (Wing ef al, 1991) and the Throughput model 


developed by the USAF/SAA (Yost, 1994). 
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In 1994 the Naval Postgraduate School conducted research in response to a 
request from the U.S. Air Force Studies and Analyses Agency (USAF/SAA) and 
culminated with the development of Throughput IJ, a model described in a 1994 NPS 
M.S. thesis by Capt. Lim Teo Weng (Lim, 1994) and enhanced the following year 
(Morton, ef al., 1995). Throughput II is a strategic airlift assets optimization problem 
formulated as a multi-period, multi-commodity flow model with a large number of side | 
constraints. It is implemented in the General Algebraic Modeling System (GAMS) 
(Brooke ef al/., 1992), and its purpose is to minimize late and non-deliveries subject to 
physical and policy constraints, such as aircraft utilization limitations and airfield 
handling capacities. 

Throughput II and other existing strategic airlift optimization models typically 
assume all data is known prior to making a decision. However, assuming that all data 1s 
known with certainty is not always realistic. Specifically, existing optimization models 
fail to properly address aircraft reliability, which is an inherently random aspect of a 
Strategic airlift system. Grounded aircraft that require repair work can significantly 
degrade the performance of an airlift system. Failing to properly model aircraft 
reliability may result in models that are too optimistic with respect to throughput 
capability. Another limitation is these optimization models may schedule unreliable 
aircraft through capacity limited airfields or airfields that have limited repair capability. 

Techniques from stochastic programming are used to develop a model that 


| incorporates aircraft reliability. Data is available for aircraft breakdown and repair rates, thus 
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allowing the development of empirical probability distributions. As a result, aircraft 
reliability can be modeled as a random variable with a known distribution. With this data, a 
stochastic extension is added to the Throughput IT model resulting in a two-stage stochastic 
linear program with recourse. The stochastic model, in addition to the features of 
Throughput II, accomplishes the following: 

e Selection of aircraft routes by anticipating potential bottlenecks in the system. 


e Minimizes the number of unreliable aircraft using capacity limited airfields or 
airfields that have limited repair capability. 


e Achieve a flow of cargo to the theater that is not interrupted by the random events 
of aircraft reliability. 


However, there is a price to be paid with respect to model size for incorporating 
aircraft reliability. Due to the large number of scenarios involved, this linear program 1s 
large in scale and beyond the capability of mieeiing languages like GAMS. To overcome 
this challenge, : special purpose algorithm designed for stochastic optimization models is 
utilized, specifically Benders decomposition. 

To objectively analyze the deployment schedule recommended by the stochastic 
program, a detailed discrete-event simulation model of the strategic airlift system 1s also 
developed. The simulation model attempts to execute deployment schedules and allows 
the user to analyze the realism of the recommended schedule. Specifically, the 
simulation model is used to execute the ‘ploiien schedule from an optimization 
solution and allows the user both to analyze the deterministic and stochastic models and 


also to compare their proposed deployment schedules. 
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l. STOCHASTIC MODELING FOR AIRLIFT MOBILITY 


Strategic airlift is a critical factor in military operations. The demand to 
distribute personnel, supplies, and equipment to combat areas can be immediate and of 
immense proportions. To project power across the world, strategic airlift, along with 
sealift, must transport the necessary units and equipment and also sustain those units. 
During Operation Desert Shield strategic airlift delivered 15 percent of the approximately 
3.5 million stons of dry cargo (Gulf War Air Power Survey, 1993, p. 6), and between 7 
August 1990 and 10 March 1991, airlift delivered approximately 23 percent of all cargo 
(Gulf War Air Power Survey, 1993, p. 112). 

This thesis examines the major issue of airlift mobility and is a continuation of 
optimization modeling research conducted at the Naval Postgraduate School in 1994. In 
that year research was performed in response to a request from the U.S. Air Force Studies 
and Analyses Agency (USAF/SAA) and culminated with the development of Throughput 
II, a model described in a 1994 NPS MLS. thesis by Capt. Lim Teo Weng (Lim, 1994) and 
enhanced the following year (Morton, ef al., 1995). This thesis describes a stochastic 
extension to Throughput II to account for aircraft reliability. A detailed discrete-event 
simulation model of the strategic airlift system is also developed so that this stochastic 
system may be more closely analyzed. This simulation model attempts to execute 
deployment schedules and allows the user to analyze the deterministic and stochastic 


models and to compare various deployment schedules. 





A. BACKGROUND 


The magnitude of the airlift effort during Desert Shield/Storm was 
unprecedented. At the height of the war during Fall 1990, the Air Force averaged 17 
million ston-miles per day of cargo and troops and by 10 March 1991, strategic airlift had 
moved more than 500,000 people and 540,000 stons of cargo (Gulf War Air Power 
Survey, 1993, p. 3). If the United States was to experience a future Operation Desert 
Shield/Storm type scenario, massive amounts of | equipment and large numbers of 
personnel would have to be transported over continents and oceans, with an impending 
deadline. The magnitude of such a deployment imposes great strains on current and 


future air, land and sea mobility systems. 


Despite the strategic success of Operation Desert Storm/Shield, the airlift 
operation experienced execution planning problems and a shortfall in airlift capability. 


Examples of each are as follows: 


e Execution. Computer models could not analyze the airlift schedule and 
determine where the flow exceeded the throughput capacity of the airfield 
structure. Consequently, too many aircraft were passing through some parts 
of the system at one time, and airfields became backlogged. As a result the 
Air Force had to halt the flow on several occasions. (Gulf War Air Power 
Survey, 1993, p. 91) 


e Shortfall. The Air Force relied heavily on the civil airline industry to fulfill its 
airlift requirements. Furthermore, without the thousands of missions flown by 
civilian air carriers, the Air Force could not have moved the required troops 
and cargo to the Arabian Peninsula by the time the United Nations deadline 
expired on 15 January 1991. (Gulf War Air Power Survey, 1993, p. 3) 





The U.S. military services are currently working to improve the airlift system, and as a 
result, various optimization and simulation models have been developed. These 
analytical tools have been used to help improve the effectiveness of limited lift assets in 
order to minimize delivery shortfalls. 

After the Gulf War in 1991, Congress commissioned a Mobility Requirements 
Study (MRS) of the United States Armed Forces. The MRS studied all aspects of 
mobility (domestic transportation, inter-theater lift, intra-theater lift and prepositioning) 
to determine the proper mix of sea, air, and amphibious lift, surface transportation and 
prepositioning. | The goal was to provide Congress with an integrated plan for procuring 
the necessary lift for power projection in the 21st century. Two linear programming (LP) 
optimizations models were developed as part of the MRS. They are the Mobility 
Optimization Model (MOM) developed by the Joint Staff's Force Structure, Resource 
and Assessment Directorate (J8) (Wing ef al., 1991) and the Throughput model 


developed by the USAF/SAA (Yost, 1994). 
1. Mobility Optimization Model 


MOM 1s a multi-period, muleccommodity model with a single source in the 
United States and a single destination in the theater of concern. This model considers 
both air and sea mobility and is used in determining the proper level and mix of lift assets 
(air, sea and prepositioning) necessary to support U.S. power projection needs. Although 
MOM serves its purpose well by modeling the spectrum of strategic lift, it is not suitable 


for answering the detailed air mobility questions currently sought by USAF/SAA. In 





short, it does not model the airlift infrastructure network and constraints in sufficient 
detail; e.g., the maximum number of planes on the ground (MOG) constraint at airfields 


is not directly modeled. 
Zi Throughput Model 


The Throughput model is a single-period, multi-commodity model which 
incorporates the strategic airlift network. This static airlift mobility model was 
developed by USAF/SAA to determine which part of the airlift system is restricting flow. 
It has also been used to evaluate alternative military strategic airlift fleets, route 
structures and basing schemes for different scenarios. Although the Throughput model 
has proven useful for USAF/SAA, it has inherent limitations due to its static nature 
(single period). For example, important constraints such as delivery time windows 
cannot be modeled when the time domain is not incorporated into the model. In 
addition, the model is unable to provide answers to key questions that concern decision- 
makers such as, on what day will Unit X commence movement and on what day will it 


be closed? 
B. THROUGHPUT Il 


Throughput II is a strategic airlift assets optimization problem formulated as a 
multi-period, multi-commodity flow model with a large number of side constraints. It 


combines the time period aspect of MOM with the infrastructure network of Throughput. 


It is implemented in the General Algebraic Modeling System (GAMS) (Brooke e/ al., 





1992), and its purpose is to minimize late and non-deliveries subject to physical and 
policy constraints, such as aircraft utilization limitations and airfield handling capacities. 
For a given fleet and infrastructure, the model can provide insight when answering many 
mobility questions, such as: 1) Are the aircraft and airfield assets adequate for the 
deployment scenario? 2) What are the impacts of shortfalls in airlift capability? 3) Where 
are the system bottlenecks and when will they become noticeable? These types of 
analyses can be useful when selecting airlift assets and investing or divesting in airfield 
infrastructure. (Morton ef al., 1995) 

In the Air Force analysis community, simulation modeling is currently used more 
widely than optimization. Simulation models more readily accommodate uncertainty and 
handle higher levels of detail, such as tracking individual airplanes and flight crews. For 
example, the Air Mobility Command’s (AMC) Mobility Analysis Support System 
(MASS), a simulation model, is a discrete-event global airlift simulation of military and 
commercial airlift assets used in strategic and theater operations. MASS is the current 
standard for all airlift studies and in recent years, no critical airlift analysis has been 
performed without comparing the results with output from MASS. Yet even MASS has 
disadvantages; it can only answer “what-if” questions, not “what’s best” questions; e.g., 
are there enough airlift assets for this scenario vs. which aircraft are the most efficient 
and versatile? Furthermore, the simulation program expends more manpower time and 
computer resources than optimization models. For instance, it requires two weeks to 


setup, run and analyze results for a typical MASS scenario. 








C. PROBLEM STATEMENT 


Throughput II and other existing strategic airlift optimization models typically 
assume all data is known prior to making a decision. However, assuming that all data is 
known with certainty is sisi always realistic. Specifically, existing optimization models 
fail to properly address aircraft reliability, which is an inherently random aspect of a 
strategic airlift system.’ Grounded aircraft that require repair work can significantly 
degrade the performance of an airlift system. Failing to properly model aircraft 
reliability may result in models that are too optimistic with respect to throughput 
capability. Another limitation is these optimization models may schedule unreliable 


aircraft through capacity limited airfields or airfields that have limited repair capability. 


D. METHODOLOGY 


Techniques from stochastic programming are used to develop a model that 
incorporates aircraft reliability. Data is available for aircraft breakdown and repair rates, thus 
allowing the development of empirical probability distributions. As a result, aircraft 
reliability can be modeled as a random variable with a known distribution. 

The stochastic optimization model is based on two-stage stochastic linear 
programming with recourse. The decision making process and uncertain events unfold in the 
following manner. First, an aircraft schedule is devised. Then, based on this schedule, the 


model examines all possible scenarios involving aircraft repair time realizations and the 


'While MASS has the capability of running with specified probability distributions for ground time, it is 
rarely run in this mode due to computational requirements (Waisanen, 1994). 





relation of these repair times to the available airfield capacities. If, for a specific scenario 
and time period the available capacity at an airfield is exceeded, a penalty in proportion to 
the excess will be incurred. The penalty function will discourage deployment schedules that 
frequently (as measured by the probability distribution) result in airfield capacity violations. 
This will produce an optimal solution that maximizes the expected value of system 
throughput—a flow of cargo from origin to destination that is not interrupted by the random — 
events of aircraft reliability. This approach represents the logic of an experienced scheduler 
in that, relative to a deterministic model, the stochastic model is encouraged to: 
e Select aircraft routes by anticipating potential bottlenecks in the system. 


e Prevent unreliable aircraft from using capacity limited airfields or airfields that 
have limited repair capability. 


e Achieve a flow of cargo to the theater that is not interrupted by the random events 
of aircraft reliability. 


However, there is a price to be paid with respect to model size for incorporating 
aircraft reliability. Due to the large number of scenarios involved, this linear program is 
~ Jarge in scale and beyond the capability of modeling languages like GAMS. To overcome 
this challenge, a special purpose algorithm designed for stochastic optimization models ts 
utilized, specifically Benders decomposition. 

To objectively analyze the deployment schedule recommended by the stochastic 
program, a detailed diceciesuenl simulation model of the strategic airlift system Is also 
developed. The simulation model attempts to execute deployment schedules and allows 


the user to analyze the recommended schedule. Specifically, the simulation model 1s 





used to execute the deployment schedule from an optimization solution and allows the 
user both to analyze the realism of the deterministic and stochastic models and also to 
compare their proposed deployment schedules. Also, examining the results of an 
optimization model via ; simulation model can provide useful feedaack with respect to the 
validity of various linear programming modeling assumptions. 

In Chapter II, the features, assumptions and limitations of the Throughput I model 
are discussed. The details of the mathematical formulation can be found in Appendix A. 
Also in Chapter II the issue of Maximum on Ground (MOG) and the impact of aircraft 
reliability on MOG is discussed. This is followed by the stochastic extension to Throughput 
II in Chapter IIT. The Benders decomposition algorithm used to solve the large stochastic 
model is described in Chapter IV. The simulation model used to compare the stochastic and 
deterministic models is developed in Chapter V and the analytical results and insight derived 


from this comparison are contained in Chapter VI. 





I. DETERMINISTIC MODEL AND AIRCRAFT RELIABILITY 


MOG (Maximum on Ground) is used extensively in the strategic airlift system 
and has a dramatic impact on throughput capability. Therefore, any detailed model of 
the airlift system must incorporate MOG. An explanation of MOG and why it 1s 
important is contained in Section A. Then in Section B, this thesis examines the model 
features and assumptions of an existing deterministic optimization model, Throughput II. 
As stated earlier, this deterministic model assumes perfectly reliable aircraft and 
deterministic ground times. But these assumptions err on the side of optimism. This 
shortcoming is also discussed in Section B. To examine the issue of aircraft reliability, 
Section C discusses the available reliability data and its impact on_ strategic 


airlift operations. 
A. MAXIMUM ON GROUND (MOG) OR AIRFIELD CAPACITY 


A key factor in any strategic airlift model is MOG. By definition MOG is “the 
highest number of aircraft being used in an operation which will be allowed on the 
ground during a given span of time based on simultaneous support.” The key phrase in 
this definition is “simultaneous support;” an airfield’s MOG value is the maximum 
number of aircraft that can be serviced at the same time when the airfield is equipped 
with specified service facilities. The MOG value for an airfield depends on a number of 


factors including parking spaces available (also called ramp space), material handling 





equipment, maintenance capabilities, fuel availability, as well as the type of service 
required (onload, refuel, or offload). Comparative scenarios are as follows: 


e Scenario A: Airfield Bravo is being used as an enroute base. In this role the 
limiting factor for airfield Bravo is the number of fuel pumps available. 
Assuming the airfield can only refuel eight aircraft simultaneously, a MOG of 
eight is established. 


e Scenario B: Airfield Bravo is an onload base. In this case the airfield may be 
limited by its ability to onload equipment due to a limited number of forklifts 
and/or cranes. Assuming the airfield can only load two aircraft 
simultaneously, this results in a MOG of two. Thus scenario B is more 
limiting than the fuel pump scenario of A. 


As noted above MOG is not static; depending on the scenario and its limiting factors, 
MOG, for example, can range from a value of eight to a value of two. 

A further variable related to MOG is ramp space. This is the maximum number 
of aircraft that the airfield can park at any given time. Thus for any airfield, its MOG 
value will be less than or equal to its ramp space. In an airlift operation, one might be 
inclined to send more planes above the MOG value because the large ramp space 
available at an airfield is believed to allow a larger flow of cargo to the theater. The 
following quote illustrates why MOG is important (versus ramp space) and why we 
should schedule in such a way to not violate MOG. 

“As more and more aircraft enter a MOG constrained system, it becomes 

progressively more difficult to serve these arriving aircraft in a scheduled 

ground time. The coordination of limited fueling, maintenance, loading 

and manpower assets to maximize throughput depends upon a limited 

number of aircraft per day. Adding aircraft can magnify the delays 

experienced for all users and throughput to a theater can and does 
decrease as the number of aircraft increases. Somalia 1s a good example. 


Early attempts to maximize flow rates into the country resulted in holding, 
diverts, malpositioned cargo, and rescheduled missions. All because too 





many aircraft were trying too operate in too constrained an environment. 

Ramp space is just the beginning. In Desert Shield, more aircraft should 

certainly deliver more to a huge ramp like the one at Dhahran: but what 

happened? Airlift aircraft waited behind commercial freighters, fighters, 
bombers, and other multinational forces waiting for limited supplies of 

fuel. Ground times increased to the point that we held missions in Europe 

for scheduled slot times down range. Slot times were not the result of 

inadequate ramp space, but were a solution to lengthy ground times 

waiting for fuel. These delays extended crew days beyond legal limits. It 

was not an effective strategy to send in a maximum flow. We achieved 

maximum throughput when scheduled slot times reduced the number of 

aircraft arriving per day.” (Merril and Szabo, 1994, p. 4) 

This quote indicates that the limiting factor in Dhahran for throughput capability is not 
ramp space but is fuel capacity (available to airlift aircraft), which defines its 
MOG value. 

As stated earlier MOG is the maximum aircraft on the ground an airfield can 
service simultaneously. But in airlift operations and in the modeling world, it 1s 
convenient at times to apply a transformation to the definition of MOG; MOG can be 
transformed to the maximum flow rate of inbound aircraft an airfield can sustain. For 
instance, if an airfield has a MOG value of two (assuming one hour ground times), this 
airfield can support aircraft arriving every 30 minutes. In military airlift operations, this 
transformation of MOG enables personnel in the Air Force to make better scheduling 
decisions. This transformation is also necessary for models with discrete time periods—in 
particular the mathematical optimization models discussed and developed in this thesis. 

Using the transformation of MOG, is it realistic to have aircraft arrive at this 
maximum flow? This question highlights two separate issues. 

1. Airlift operations. The current planning factors require an aircraft be serviced 


in an allotted amount of time for refueling, offloading or onloading. But what 
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happens when an aircraft breaks down? To remedy the situation, additional 
manpower or equipment resources will be required. But from where? Since 
the airfield has no additional resources readily available, the only solution is 
to divert resources from other aircraft being serviced, thus multiplying the 
problem situation and decreasing capacity even further. And as more aircraft 
arrive, significant delays can develop while the base is operating beyond its 
current capacity. 


2. Modeling. In building a mathematical model with discrete time steps, MOG 
is transformed from a physical capacity (i.e., 2 plane slots) into a unit that 
involves time (i.e., 48 plane slots per time period; assuming a 24 hour time 
period and one hour ground times). From the perspective of the model, it is 
perfectly acceptable to have all aircraft (i.e., 48 aircraft) arrive at one time. 
However, to fully utilize this capacity in reality, the model must produce a 
perfect back-to-back schedule (i.e., aircraft must arrive every 30 minutes). If 
aircraft arrivals are not uniformly spaced, even if aircraft do not break down, 
the model cannot achieve this capacity in reality, even though the model 
thinks it can. 

To remedy this potential problem; one should lower the arrival rate or flow of 
aircraft to this airfield from the maximum. This will result in 1) more manpower and 
equipment immediately available for broken aircraft and 2) minimize the possible 
violation of MOG in the optimization model due to discrete time periods. This will 
allow aircraft to land and immediately be serviced, thus achieving a realistic flow of 
cargo from origin to destination. This thesis will determine, by modeling one major 
aspect of ground time, namely aircraft reliability, the proper rate of flow of aircraft from 


origin to destination. This flow rate can then be transformed to determine at what 


percentage of MOG an airfield should operate. 
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B. THROUGHPUT II 


The original Throughput II model was developed by Capt. Lim Teo Weng (Lim, 
1994). Since then the model has been enhanced in ongoing research at the Naval 
Postgraduate School (Morton, et al., 1995). A brief overview of the mathematical 


formulation of Throughput II is in Appendix A. 
1. Model Features 


The Throughput II model is a multi-period model for determining the maximum 
on-time throughput of cargo and passengers that can be transported with a given fleet 
over a given network, subject to appropriate physical and policy constraints. The 
objective function minimizes the total weighted penalties for late deliveries and non- 
deliveries. This model considers only inter-theater, not intra-theater deliveries. The 
major features -of the airlift system currently captured by the model include (Morton, ef 
al, 1995), 

e Multiple origins and destinations. The model flies aircraft through multiple 
origin, enroute and destination airfields. Due to the size of the airlift system, 
ageregation of origin and destination airfields is done to allow tractability of 
the model. | 

e Flexible routing structure. The air route structure supported by the model 
includes delivery and recovery routes with a variable number of enroute stops 
(usually between zero and three). This will give the model the option of 
short-range flights with heavier loads or long-range flights with lighter loads. 


For further routing flexibility, the model also allows the same aircraft to fly 
different delivery and recovery routes on the same mission. 


e Aircraft-to-route restrictions. The user may impose aircraft-to-route 
restrictions; e.g., military aircraft may only use military airfields for enroute 
stops. This particular provision arises because the USAF AMC may call upon 








civilian commercial airliners to augment USAF aircraft in a deployment, 
under the Civil Reserve Airfleet (CRAF) program. The model distinguishes 
between USAF and CRAF aircraft. 


Aircraft assets can be added over time. This adds realism to the model, 
because CRAF and other aircraft may take time to mobilize and are typically 
unavailable at the start of a deployment. 


Delivery time windows. In a deployment, a unit is ready to move on its 
available-to-load-date (ALD) and has to arrive in the theater by its required- 
delivery-date (RDD). This aspect of the problem has been incorporated in the 
model through time windows that are specified for each unit. The model 
treats the time window as “elastic” in that the cargo may be delivered late, 


subject to a penalty. 


Assumptions 


The assumptions used in the model are as follows (Morton, ef al., 1995): 


Congestion. Inventoried aircraft at origin and destination airfields are 
considered not to affect the aircraft handling capacity of the airfield. This 
assumption is not strictly valid since an inventoried aircraft may interfere with 
the ability of the airfield to process aircraft in the airlift operation. 


Deterministic ground time. Aircraft turnaround times for onloading and 
offloading cargo and enroute refueling are assumed to be known constants, 
although they are naturally stochastic. This ignores the fact that deviations 
from the given service time can cause congestion on the ground. To offset the 
optimism of this assumption, Throughput II uses an efficiency factor (<1) in 
the formulation of MOG consumption constraints to lessen the impact 
of randomness. 


Computational tractability. The model aggregates onload and offload 
airfields, and uses continuous decision variables and discrete time periods— 
typically one day. 








Limitations 


Throughput II errs on the side of optimism with respect to throughput capability 


due to the assumption of deterministic ground times and perfectly reliable aircraft. Also 


the single-day time periods can decrease the model’s fidelity. This results in the 


following potential problems that will be examined in this thesis: 


Deterministic ground times for aircraft. The model, if needed, will 
continually saturate airfield capacities to achieve the maximum flow of cargo 
from origin to destination. However, the aircraft ground times used in the 
calculation of MOG consumption represent the expected times for 
onload/refuel/offload, resulting in optimistic throughput capability (via 
Jensen’s inequality). In an attempt to overcome this optimism, the 
deterministic model uses a MOG efficiency factor to lessen the impact of 
random ground times. 


Perfectly reliable aircraft. Broken aircraft can have an immediate impact on 
throughput capability, especially when airfield capacity resources are limited. 
Moreover, the necessary repair work consumes resources that might be 
needed elsewhere. The location of the broken aircraft is also critical. If an 
aircraft breaks enroute, not only does the necessary repair work consume 
resources at the enroute airfield, the cargo onboard the aircraft will be delayed 
in arriving at its destination. If the aircraft breaks at an airfield with limited 
maintenance capability, the required maintenance facilities and assets might 
not be readily available, resulting in even longer delays. 


Time resolution. With one-day time periods, the model can route aircraft in a 
manner that causes local congestion. For example, all aircraft could arrive at 
an airfield within a small time window instead of being dispersed over an 
entire day. In reality this would cause local congestion, even though the 
model’s representation of aircraft handling capacity is observed. Another 
limitation is that the model rounds the time (to the nearest day) at which an 
aircraft arrives at origin and destination airfields. 
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C. AIRCRAFT RELIABILITY 


There are many factors that contribute to the stochastic nature of ground times for 
aircraft. Some of these factors include local congestion at an airfield, lack of 
maintenance personnel or spare parts, airlift users, aircraft reliability, and weather. 


Aircraft reliability accounts for a significant portion of random ground times. 
1. Data 


The Air Mobility Command (AMC) collects aircraft breakdown and repair rates 
for each of the strategic military aircraft, which are defined as follows (AMC Pamphlet 
21-2, p. 17): 


e Break: system malfunction that renders aircraft non-mission capable (NMC) 
after landing. 


e Break rate: percent of aircraft landings that have system malfunctions 
rendering aircraft NMC. 


e Break rate formula: 


BUCS nt (2.1) 
number of landings 


e Fix: completing maintenance actions on NMC aircraft rendering the aircraft 
airworthy (PMC, partially mission capable, or FMC, fully mission capable). 


e Fix rate: percent of aircraft landing NMC that are fixed within established 
time frames ( 0-4, 4-8, 8-12, 12-16, 16-24, 24-48, 48-72, >72 hours). 


e Fix rate formula: 











number of aircraft fixed within a specified time frame 
100. (2.2) 
number of aircraft landing breaks 


Appendix B contains data in this form from 1994 for the C5, C17 and C141. 
Based on this data, an empirical probability mass function, p,(4), for each aircraft type a 
is developed. Each possible realization of the random variable 7,, representing the 
required repair time, will have an associated probability, P(7, = 4 = p,(t). This 
probability mass function applies each time a plane lands to determine the required repair 
time, 7,. This assumes that the probability mass function does not depend on the hours 
flown by the aircraft (This assumption is not entirely correct, since aircraft reliability 
typically decreases with increasing flying hours.). With this assumption, the probability 


mass functions are derived resulting in nine realizations per aircraft (Table U-1). 


— Pi Dcs(t) Dc17(t) Do141(t) 


0.86540 0.93110 
0.03078 0.02733 





0.02933 0.00832 
0.01643 0.00238 0.02116 
0.01582 0.00238 0.01205 


0.01269 | 0.00832 0.00910 
0.01872 0.01069 








Table I-1: Probability Mass Function, p,(/). 
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In review, we have the following notation: 
e 7,,-tandom variable representing the repair time of an aircraft of type a. 
¢ p(t) - the probability mass function for 7,,. 


e S.- support of the random variable 7,,, where S, = {0, 2, 6, 10, 14, 20, 36, 60, 
72. 


To consider aircraft reliability, a model can be formulated that tracks each 
individual aircraft. By using binary variables, the model can determine if that specific 
aircraft breaks when it lands at airfield 5 in time period t, thus the binary variable will 
require three indices (aircraft identification number, airfield and time period). But 
because the model does not know a priori how to schedule the aircraft, the model does 
not know which airfield the plane will land at on any given time period. Therefore the 
model must consider every possible scenario that can occur, Fora singe aircraft, the size 


of the sample space (assuming nine realizations per aircraft) is 


| = Q (total number of airfields)*(time periods) (2.3) 


In extending this to all aircraft in the model, the size of the sample space (assuming 9 


realizations for each aircraft) is 


IQ = 9 (total number of aircraft)*(total number of airtields)*(time periods) (2.4) 
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To examine these scenarios in a mathematical program, the model will require a 
constraint for each scenario (as a minimum); and a binary variable for each individual 
aircraft (as an indicator variable) to determine if that particular plane lands at airfield 5 in 
time period ¢. Thus the size of the model is proportional to the size of the sample space. 


Using the Throughput II data considered by Lim (1994), this results in a minimum of 


; ‘ eee 
| Ql = g (416 aircraft)*(16 airfields)*(30 time periods) == g 199,680 (2.5) 


constraints, in addition to the requirement for binary variables. A model involving this 


minimum number of constraints is impossible to solve given today’s technology. 


Z: Scenario Development 


To reduce the size of the sample space and the associated model, three steps are 
enacted. First, the convolution is used. For instance, the convolution of p,(t) with 


itself is 


Palt)= > p,(t- y)p,(y). (2.6) 


yes, 


The convolution of p,(f) with itself n times is denoted p,,(f): the probability mass 


function of the total repair time random variable, 7,,,, of n aircraft of type a. Inthe model 
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we develop, the convolution procedure will be utilized to reduce the dimension of the 
model. In review, we have the following additional notation: 


e 7, - random variable formed by summing n independent and identically 
distributed random variables, 7,, each with probability mass function p,(‘). 


® Pat) - the probability mass function for 7, with support S,,. 


To reduce the number of scenarios even further, the second step approximates the 
distribution of he random variable, 7, with another random variable, ie . This new 
random variable consists of a smaller number of realizations. While T., has fewer 
realizations than 7,,,, the probability mass function p, (rt) still closely resembles that of 


Pat). This results in the following notation: 


A 


e 7, ~tandom variable that approximates 7, but has fewer realizations. 
e p(t) - the probability mass function for i with support Ss. 
The third step incorporates the random variable ie in the stochastic model 


developed in this thesis. We define the random variable 
RTime, = +7; (2.7) 
Hi 


realizations of this random variable are denoted RTime®. Multiplication of this random 


variable by the number of landings of aircraft of type a gives the total repair time 


required; this has the correct distribution if exactly n aircraft have landed. It will be 
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shown in Chapter III that the stochastic model developed examines separately the 
possible scenarios for each base and time period, 1.e., each (5,4) combination will have its 
own random variable RTime,. Thus the correct notation is RTime,,,. In review we have 


the following notation: 


e RTime,,, - the repair time random variable for aircraft of type a at airfield 5 
and time period f. 


¢ p,,,(t) - the probability mass function for R7ime,,, with support ee 
These three steps reduce the size of the sample space for each base and time 


period combination, Q,,, to (assuming nine scenarios per aircraft type (7 = 9) for p,,,(¢) ) 


Q,,| = 9° aircraft types) (2.8) 


This results in a minimum of an additional 
g > arcralt types) x (16 airfields)*(30 time periods) = 349,920 (2.9) 


constraints. The stochastic model developed in Chapter IJI will then require a minimum 
of 349,920 constraints. The advantages and disadvantages of these three steps are 
as follows: 


Advantages: 
e Reduces the number of scenarios or the size of the sample space. 
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No need for integer programming. Since R7ime,,, approximates the repair 
time variable for a single aircraft of type a, at airfield 5, time period f, 
multiplication of R7ime‘, by a continuous variable gives the approximate 
total repair time required in scenario @ (assuming the effects of using a 
continuous decision variable are negligible). 


A stochastic model using this sample space is too large to be solved using 
modeling languages like GAMS, but it can be solved using decomposition 
based algorithms. 


Disadvantages: 


Since the input to the stochastic model will be R7Time,,,, the realizations of 


this random variable will be exact if the model sends exactly n planes of type 
a to airfield 5 in time period ¢. In this case, the only error is caused by using 


A 


7, im approximating 7,,. This error will get smaller as the number of 


an 


realizations per aircraft type used in developing ji approaches the number of 


realizations of 7,,,. 


If the model sends more or less than n planes, the realizations of i is not 
accurate since this was based on nv planes. This error is in addition to the error 


caused by using T., to approximate 7,,,,. 


D. EFFECT OF AIRCRAFT RELIABILITY ON AIRLIFT OPERATIONS 


In developing the stochastic model, two separate modeling issues were of 


concern. First, a broken aircraft consumes a portion of ramp space and second, a broken 


aircraft will have some type of impact on MOG. Consider the following points relating 


to the issue of broken aircraft consuming ramp space: 


In the literature ramp space is never mentioned as the limiting factor for 
throughput capability. During the Gulf War, Rhein-Main had sixty-eight 
aircraft on the ramp, even though the base had parking spaces for only fifty- 
six aircraft (Gulf War Air Power Survey, 1993, p. 100). 
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e Ramp space data is not currently available to support a stochastic model with 
this model enhancement. 


In light of these two points, this thesis does not attempt to model the impact of 
broken aircraft on ramp space, i.e., the stochastic model does not have a ramp space 
consumption constraint. On the other hand, consider the following points: 

e A broken aircraft does consume ramp space, but the plane also requires 
resources to complete repairs—resources that could be used for processing 
planes in the delivery of cargo. 

e In peacetime there is no urgent need to repair aircraft. But during a critical 
airlift operation, a broken plane has an immediate impact when aircraft supply 
is a limiting factor, the aircraft is carrying a critical load or the aircraft breaks 
at an airfield with limited maintenance capability. In these cases additional 


resources may be required to repair the plane. 


e In an airlift operation the goal is to achieve the highest throughput capability 
possible, but this is directly related to MOG and not ramp space. 


Thus, the stochastic model must mathematically model the impact of aircraft reliability 
on MOG. | 

Using another transformation of MOG, this thesis will be referring to hours of 
MOG consumed. For instance, if an airfield has a MOG of two, then this airfield has 48 
hours of MOG available (24 hours per day * 2) per day. So when a plane does break, 
how many hours of MOG does the plane actually consume? In addition to the resources 
required for the onload/refuel/offload, an additional amount of resources is required to 
repair the aircraft. Converting repair times directly into MOG consumption is often too 
pessimistic; instead only a fraction of the total repair time actually reduces capacity. The 


longer the aircraft is on the ground the smaller that fraction. If the repair time random 
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variable realization is two hours, for example, we assume that all repairs will consume 
an additional two hours of MOG. On the other hand, if the repair time random variable 


realization is 72 hours, we assume the following: 


e If a plane breaks at a large enroute airfield with unlimited maintenance 
personnel, an additional six hours of MOG are consumed. This assumes six 
hours are required to offload and then in turn onload the cargo onto another 
plane, while the original plane is repaired at the hangar. 


e If the plane broke at an airfield with very limited maintenance personnel, the 
plane will consume seventy-two hours of MOG. 


Using these ideas, the following is performed to produce data used in the stochastic 
optimization model for each base and time period combination: 


1. Obtain the probability mass function, p,(t), for the random variable 7, for 
each aircraft of type a. 


2. Modify the realizations of 7, to account for the amount of MOG hours 
consumed. This can depend upon the airfield 5 and the time period ¢. 


3. Perform the convolution and approximation to obtain the random 
variable C. 


4. Divide the realizations of 7, by n to obtain the random variable MRT" Ime ,, . 
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lil. STOCHASTIC OPTIMIZATION MODEL 


Due to the variety of stochastic optimization models discussed in the literature, 
Section A gives a broad overview of stochastic programming terminology related to the 
stochastic optimization model developed in this thesis. Then in Section B, the stochastic 
extension to Throughput II is described, resulting in a two-stage stochastic linear program 
with recourse. This will permit a study of how the optimal solution changes when going 
from a deterministic model to a stochastic model that incorporates aircraft reliability 
(Chapter VI). It 1s shown in Section C that this stochastic optimization model is very 
large and must be solved using a decomposition based algorithm. The chapter then 
concludes with a summary of the underlying constructs used to develop the stochastic 
model and the additional model features gained (in addition to those of the Throughput II 


model) by using this stochastic extension. 
A. STOCHASTIC PROGRAMMING 


In deterministic linear programs (LP), all parameters are treated as known 
constants. This can lead to an optimistic or pessimistic objective function value, 
depending on which parameters are random. Moreover, in general it may lead to sub- 
optimal decisions because they fail to hedge against the full range of possible scenarios. 
On the other hand, stochastic programming models treat a subset of the parameters as 


random variables with a known distribution. Even when the “exact” distribution of these 
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random variables can not be obtained, the stochastic model will be an improvement over 
the one-point forecast of the deterministic LP. 


In the stochastic programming models presented here, realizations of random 


parameters are denoted with superscript @, the corresponding sample space by Q, and the 
probability mass function by p". There are a number of options available in extending 


the following linear program 


minimize cx: (3.1) 
subject to Ax =b 


x>0 


to the stochastic domain. One such option is to insist that a feasible decision be made no 


matter what scenario might unfold: 


minimize Cx : (3.2) 
subject to A* x > b° VoeQ 
x20. 


Here we have assumed that the objective function coefficients, c, are known with 
certainty. Many decision makers regard (3.2) as too pessimistic a formulation because it 


1s typically driven by the most “extreme” scenario. In most real-world systems, many of 
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the structural constraints are elastic and may be violated at a certain cost. There are a 
variety of stochastic programming models that differ primarily in how this issue of 
“infeasibility” is handled. Those relating to the stochastic optimization model developed 


in this thesis are described below (Morton, 1994). 
1. Simple Recourse Models 


In this type of model, infeasibility is accepted, but only at a certain cost; with 
each constraint violation, there is an associated cost. Let g = (¢1,¢2,93 ---.¢m) denote the 
unit cost of violating each constraint. The positive part of a scalar variable, s, is defined 
as s' = max{s,0}; the positive-part-operator may be applied to a vector, componentwise. 


With this definition, the simple recourse model 1s 


minimize,,5 cx+ EF, h(x,@) (3.3) 


where h(x,@)=q(b° — A®x)’, 


where E,h(x,@) = ¥ p°h(x,@). 


wmeQ 
Zz: Recourse Models 


In recourse models infeasibility is corrected, afterwards, at a certain cost. The 


costs due to the amount of violation in the constraints are determined after the 
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observation of the random data and are denoted as recourse costs; the expected value 


criterion is used. In this case the model is written as 


minimize cox+ FE h(x,@) (3.4) 
x20 
where K(x,®) = minimize fv 


subject to Wy 2=b® — A®x 


As in the simple recourse case, E.,,h(x,@) is called the recourse function; W is called the 


recourse matrix. The recourse variable, y, measures the corresponding violation in the 
constraints, if any. If W = J, we have the simple recourse model. 
The “standard” two-stage recourse model replaces the first stage constraints x > 0 


with general structural inequalities and otherwise generalizes the problem as follows: 


minimize {cox + E h(x,@)}. (3.5) 


x20 


The second-stage cost we will incur, if we observe scenario @ and act optimally given 


first-stage decision x, is /(x,@) and may be expressed as 
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h(x,@) = minimize fey | (3.6) 


subject to W°y=d° + Bx 


y= 0. 


The recourse aspect of the model gives it a unique characteristic. The model is. 
two-stage: first the decison variables are chosen, then the stochastic variables are 
observed. These two steps will determine the recourse variables (to recover feasibility) — 
and their corresponding penalties. These penalties depend on both the constraint violated 
and the magnitude of violation. The recourse cost has the effect of restraining 
“ageressive” choices of decision variables if the costs involved with regaining feasibility 
outweigh the benefits. Thus, the two-stage recourse model minimizes the sum of our 


original first-stage costs and the expected recourse costs. 
B. MODEL 


As stated in Chapter IJ, the deterministic Throughput IT model minimizes late and 
undelivered cargo subject to system constraints. One set of these constraints, MOG 
consumption (see Appendix A, constraint A.15), contain inherently random coefficients 
because the ground times are of a random nature. This 1s where stochastic programming 
can be used. If the pessimistic formulation (3.2) is used, it would be driven by the largest 
conceivable ground time; that way we could ensure that we will never exceed MOG. 


However, this is an extreme and unacceptable solution because the airlift system would 
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be dramatically underutilized. Thus instead we will permit certain scenarios to violate 
the MOG constraint, but add an objective function formulation that discourages the 
model from generating schedules that yield frequent or large violations (3.6). 

In this section we describe a stochastic extension to Throughput II to include 
aircraft reliability. In the following pages, the first-stage variable x is a large vector that 
includes all decision variables of Throughput II (1.e., Xyars Yor, etc.). The following is 


added to the Throughput II model, resulting in a two-stage stochastic linear program model 


with recourse. 


1. Additional Data 


MOGPEN,, This penalty is used in the objective function when airfield 
(base) b, in time period ¢, exceeds MOG, referred as MOGCap;; 
these correspond to the objective coefficients, f, in (3.6). 


MRTime*,, The modified repair time random variable for aircraft type a, at 
airfield b, at time period ¢, is denoted MRTime,,,. Realizations of 
this random variable are denoted MRTime*,. These values 
correspond to coefficients in the B® matrix in (3.6). 


2. Additional Decision Variable 

R, The amount by which MOGCap, 1s exceeded at airfield 6 in time 
period ¢ for scenario @; this variable corresponds to the y 
variable in (3.6). 


3. Additional Objective Function Term 


The sample space for the vector of the modified repair times for all aircraft types 


at airfield } at time period 1s denoted Q,,={1, 2, ..., Ks,}. From that sample space, each 
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outcome w(b,t)<¢Q,, has an associated probability p;)°. This results in a recourse 


function for each (6,f) combination with the following notation 


Ee o(b,t vt) Pon (%, o(b,t)) = DP Pal o(b,t)), (3.7) 
ab, 1)eQs, 
where 
h, (x, @(b,t)) = MOGPEN,,* R,, - (3.8) 


In the following pages this thesis abbreviates w/b,t) by @, resulting in the following 


Eghy(x,0) = >) Ply, (2,0). (3.9) 


oeQ,, 


Using this notation, the following term is added to the original objective function of 


Throughput II: 


YY Euhy(2,.0)=¥¥ ¥ pp(MOGPEN,,* Re). (3.10) 


b i we, 
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This will penalize the model whenever the MOGCap, at airfield d, in time period 


f, scenario @ is exceeded. How the value of MOGPEN,, is determined and its impact on 


the model is discussed in Section D of this chapter. 
4, Modification to Throughput IT Constraints 


Throughput I uses the deterministic ground time for calculating the time required 
for an aircraft to execute a route and for the amount of MOG consumed. For example, to 
calculate the time a loaded aircraft will arrive at the second enroute airfield, the model 
sums the onload ground time, one enroute ground time and the required flight time. 
Using the deterministic ground time is an optimistic assumption. Aircraft reliability, 
only one aspect of the ground time variability, can cause delays in the delivery of cargo. 
Since it is unreasonable (from a tractability perspective) to have aircraft actually “break” 
in a linear program, the model needs to indirectly model this issue. This can be achieved, 
in part, by using the expected value of the modified repair time, in addition to the 
deterministic ground time, for the calculation of MOG consumption and the time for an 
aircraft to execute a route. Even though the stochastic model uses the expected value of 
the modified repair time in the MOG consumption constraint, the optimal solution is 
optimistic (via Jensen’s inequality). This optimistic solution is “modified” by the 


addition of new constraints. 
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5, Additional Constraints 


In calculating the MOG consumption at airfield 5 in time period ¢, the Throughput 
qT model uses deterministic ground dines (see Appendix A, constraint A.15). As stated 
earlier, this will lead to an optimistic solution. To account for a portion of the random 
ground times, aircraft reliability is incorporated into this constraint by using the random 
variable MR7ime,,,. The original constraints are maintained (for reasons to be discussed 


in Chapter IV) and the following constraints are added: 


1) YY Yl MoGReq,,*(GTime,, + MRTime’s, ) | 24} X arp (3.11) 
u oa reR, Velie 
1’+| DTime,,, |= 


+> Yl MocReg,,*(GTime,, + MRTime’, ) /24]*Y, 1 


a reR, rer. 
t’+| Dimes, |=t 


—~R° < MOGCap, V to 
(2) R,, 20 V b,t@. (3.12) 
The stochastic model will examine each scenario @, at airfield 6, time period ¢, and 


determine if the MOGCap, has been exceeded. If it has, R,, will be the amount of 


the violation. 
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C. AN EQUIVALENT REPRESENTATION OF THE MODEL 


The deterministic formulation of Throughput IJ is a linear program and hence can 


be written in the following standard form: 
minimize cx (3.13) 
subject to Ax =b 

The stochastic extension to Throughput II results in the following model: 


minimize at YY Eh, (2,0) (3.14) 
b ft 


subject to Ax =b 


x20, 


where h,.(x,@) = minimize MOGPEN,,* Ry, 


subject to 
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YY DY [MoGReq,,*(GTime,, + MRTime’,) / 24}* X yy 
Vel 


ua reéR, 
t’+| DTime,,, |=t 


+¥E  ¥|MoGReg,,*(Grime,, + MRTime®, )/ 24} Yup 


r rel 


— R° < MOGCap, 
(2) R, 20. 

Upon close examination of the stochastic extension to Throughput II, the 
additional constraints only examine the effect of the modified repair time random 
variable at airfield 6 and time period ¢. Thus the sample space for all aircraft at airfield b 
in time period /, Q,,, has a size of (see Chapter II, Section C) 

2,,, = 9° smn or) = 9° = 729, (3.15) 
resulting in an additional! 


729* (16 airfields)* (30 time periods) = 349,920 (3.16) 


constraints. Combined with the Throughput IJ data considered by Lim (1994), we have a 


total of 
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(349,920 + 7,671) = 357,591 | (3.17) 


constraints. The stochastic model also has 362,166 decision variables. This results in a 
very large model which cannot be solved using standard modeling languages like GAMS. 
However, it has a special structure that can be exploited by decomposition based 
algorithms, specifically Benders decomposition (Chapter IV). Using this algorithm, the 
stochastic optimization model is solved under 20 minutes on a IBM RS6000 model 590 
workstation (the deterministic model, using GAMS, solves under two minutes), obtaining 
a near-optimal solution. 

With the solution from the stochastic optimization model, this thesis examines the 
impact of aircraft reliability on strategic airlift. Specifically, the thesis examines the 
change in the optimal solution when going from a deterministic model to a stochastic 


model that incorporates aircraft reliability. The analysis is contained in Chapter VI. 


D. CALCULATION OF PENALTY 


The number used for MOGPEN,, can have a dramatic impact on the optimal 
solution. If MOGPEN,, = 0 (for all airfields b and time periods £), the stochastic model 
will give the same solution as Throughput II (ignoring the use of E[MRTime,»] ). At the 
other extreme, setting MOGPEN,, = +00 (for all airfields 6 and time periods #) will not 
allow the model to exceed MOGCapy, in any scenario; i.e., we recover the pessimistic 
model (3.2). For example, if the worst case scenario had all aircraft break for 22 hours 


(assuming a 2 hour onload per aircraft and MMOGCap, = 2), then the stochastic model will 


—_ 
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only send two aircraft to the airfield. This illustrates that the value used for MOGPEN,, 
has a dramatic impact bs the optimal solution. 

As stated in Chapter II, one goal of the stochastic model is to determine the 
optimal rate of flow of aircraft from origin to destination, a flow of aircraft that is not 
interrupted by the random events of aircraft reliability. If the stochastic model exceeds 
MOG for a specific scenario, the model needs to incur a cost for attempting to use more 
resources than available. And this cost needs to be large enough so that the model 
chooses the penalty for not delivering the cargo versus the penalty of exceeding MOG. 
Specifically, when the model sends aircraft a to airfield 5, this aircraft will require the 
following amount of resources to service the aircraft (assuming the aircraft does 


not break): 


resources needed = MOG Req,,*GTime,,,. (3.18) 


Now if airfield 6 has no additional resources available for this aircraft (i.e. it is operating 


at or above MOGCap,), then the recourse variable R,’ will increase by 


R,, = MOGRegq,,*GTime,,,. (3.19) 
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The objective function will now increase by 
MOGPEN,,* Ry, = MOGPEN,,* MOG Req,,*GTime ,, . (3.20) 


This increase tn the objective function must be larger than the maximum penalty that will 


be incurred for not delivering the cargo, 


Maximum penalty incurred for not 


MOGPEN,,* MOG Req,,*GTime,,, > I (3.21) 


delivering one plane - load of cargo 


The above equation can now be solved for MOGPEN,,. 
E. SUMMARY OF UNDERLYING CONSTRUCTS 


In developing the stochastic optimization model, various techniques were used to 
reduce the number of scenarios and to model the impact of aircraft reliability on strategic 


airlift operations. 


e Modified repair time - impact of aircraft repair maintenance on MOG. The 
actual impact 1s not accurately known and will depend on the experience of 
Air Force personnel. 
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Convolution. This technique is used to reduce the number of scenarios in the 
model. To determine n, the optimal solution from Throughput II 1s used as a 
starting point. 


ae The accuracy lost in using [- to approximate 7), will be examined. As 


the number of realizations in 7, increases, the accuracy of the approximation 
increases, but so does the number of constraints 1n the stochastic model. 


MOGPEN,,. This penalty indirectly reduces the MOG value for an airfield 


by forcing the airfield to operate at a percentage below MOG, allowing the 
airfield to have resources available to repair broken aircraft. 


E[Modifited Repair Time]. In an airlift operation there will be delays in the 
delivery of cargo caused by the unexpected grounding of aircraft for repairs. 
Using the expected modified repair time in the calculation of determining the 
arrival time of aircraft should improve the model. 


There is no method available to check the first construct. On the other hand, the 


remaining constructs can be verified, not individually, but as a whole by a simulation 


model. This simulation model will bé discussed in Chapter V. 


F. MODEL FEATURES 


This stochastic model, in addition to the features of Throughput II, accomplishes 


the following: 


Selection of aircraft routes by anticipating potential bottlenecks in the system. 


Minimizes the number of unreliable aircraft that use capacity limited airfields 
or airfields that have limited repair capability. 


With the proper choice of MOGPEN,,, the optimal solution will achieve a 
flow of cargo from origin to destination that is not interrupted by the random 
events of aircraft reliability by forcing an airfield to operate at a percentage 
of MOG. 
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iV. BENDERS DECOMPOSITION 


As shown in Chapter III, the stochastic optimization model is very large. To 
solve this optimization model, we use Benders decomposition to convert the large 
problem into many appropriately coordinated subproblems of manageable size. Section 
A gives a broad overview of how Benders decomposition works, followed by the theory 
(Van Slyke and Wets, 1969) in Section B. To solve the stochastic optimization model, an 
algorithm based on Benders decomposition theory is implemented. An explanation of 
this algorithm is contained in Section C. Then in Section D, the performance of this 


algorithm in solving the stochastic optimization model is examined. 
A. INTRODUCTION 


The stochastic optimization model is inconvenient to solve directly due to the 
large number of decision variables and structural constraints. Hence, one approach is to 
adopt a ciation Strategy in which only a few of the constraints are explicitly 
maintained. Using this approach, the strategy of the decomposition procedure is as 
follows. The large LP 1s initially separated into two LPs which are called the Master and 
the Subproblem. The Subproblem, in turn, separates nto one subproblem for each 
scenario. The Master consists of the general constraints (the deterministic constraints) 
and the scenario-subproblems consist of special constraints (the stochastic constraints). 


The decomposition algorithm begins by solving the Master first; the Master then passes 
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down a new set of right-hand-side resources to the Subproblems; the Subproblems are 
solved and then pass the Master a new set of constraints derived from these resources. 
The Master is solved again (now consisting of the deterministic constraints and a set of 
constraints derived from the stochastic subproblems) and the process repeats itself until 


an optimal solution is obtained. This process is also called a row generation technique. 


B. THEORY 


To illustrate the ideas underlying Benders decomposition this section considers 
the special case of a two-stage stochastic linear program with recourse in which there is 
only one scenario. The extension to the more general case of many scenarios 1s 
straightforward. We will also assume the second-stage subproblem is always feasible 
given any first-stage decisions (called complete recourse). The assumption of complete 
recourse 1s realistic whenever second-stage infeasibilities can be modeled via penalty 
costs. To apply Benders decomposition, complete recourse is not required, but will 


simplify the discussion. With this special case, the optimization model of interest is 


minimize cx + fy (4.1) 
subject to Ax =) 

—Bx+Wy =d 

x20, y20 
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This problem (4.1) is equivalent to 


minimize cx + h(x) (4.2) 


subject to Ax =b 


where h(x)= minimize fy (4.3) 


subject to. Wy=d+ Bx 


By taking the dual, the recourse function (4.3) can be written as 


h(x) = maximize (d+ Bx) (4.4) 


subjectto aAWs<f. 


Another assumption used in Benders decomposition theory is that the second-stage 
problem (4.3) is bounded. If it were unbounded for some finite first-stage decision, then 
by duality theory the corresponding dual must be infeasible. Since the dual feasible 
region is independent of x, unboundedness can occur in (4.3) only if it occurs for all x. 


Since the second-stage problem is bounded, the recourse function attains an optimal 
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solution at an extreme point. By enumerating the extreme points of the dual feasible 


) 


region as z") 2... 70°”, the recourse function (4.4) is equivalent to 


h(x) = maximize (d+ Bx). (4.5) 


Upon closer examination of the recourse function (4.5), it can be shown that for each — 


extreme point 2” of the dual feasible region, 2°d+2Bx, as a function of x, is a 


hyperplane; 2d is the intercept and xB is the gradient. Each one of the Z pieces of 
h(x) is called a cut and these cuts are supports (or “pieces”) of the convex piecewise 
linear recourse function. In summary, the recourse function is a piecewise linear, convex 
function generated by the extreme points of the dual feasible region. 

As a result of (4.2)/(4.5), the original linear program (4.1) can be written in a 


form that suggests Benders algorithm: 


minimize cx +0 | (4.6) 
subject to Ax =b 


cuts -G'x+0@ 22’ i=1,2,...,0 


where G' =7'"B and g! =2'd. 
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As stated earlier, Benders decomposition is an iterative process. Initially, the 
Master consists of the deterministic constraints and no cuts while the Subproblem is the 
recourse function. In the first iteration, the Master problem is solved and passes a new 
set of resources to the Subproblem. For the given first-stage decision vector, x, the 
Subproblem will pass up a new cut to the Master problem which represents a “piece” of 
the piecewise linear function A(x) that is a support at x. Now the Master problem consists - 
of the original deterministic constraints plus this cut. The process repeats itself until 
optimality is achieved. The algorithm will be computationally efficient when only a 
small subset of the potentially very large number of cuts needs to be generated. 

During any iteration of Benders algorithm, only a subset of the necessary cuts 


have been added to the Master problem. At each iteration, the Master problem is solved 


for optimality, denoted (%,6) ,where the scalar value @ is an estimate of the recourse 


function, A(x). If the set of cuts in tie Master problem is complete at < then @ = A(Z) 
and the algorithm terminates with an optimal solution. On the other hand, if the set of 
cuts in the Master problem is incomplete at x, then @< h(x) and a cut is added to the 
Master problem. As the set of cuts in the Master problem is augmented each iteration, 
“more informed” first-stage decisions are made. The algorithm terminates when the 
Master has sufficient information regarding the recourse function in the neighborhood of 
an optimal first-stage decision. 

During the execution of Benders algorithm, an upper and lower bound on the 


optima! objective function value can be calculated. At any given iteration, the set of cuts 
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in the Master problem provide a outer linearization of the convex recourse function. 
Therefore, the lower a on the objective function value is given by cf + 7 As each 
iteration progresses, the Master problem becomes more constrained as cuts are added, 
therefore the lower bounds generated in successive iterations are increasing. On the 
other hand, the upper bound on the objective function value is given by cf+ ff. This 
upper bound will not necessarily decrease in successive iterations, therefore the 
algorithm should save the solution that corresponds to the least upper bound generated 


so far. 

At any given iteration, if the lower and upper bounds of the objective function 
value are equal, the algorithm terminates. On the other hand, if the upper and lower 
bounds do not match, then the Subproblem passes a new cut to the Master problem and 


another iteration is completed. 
C. ALGORITHM 


It was shown in Chapter III, Section C that the stochastic optimization model 


developed in this thesis has the following form: 


minimize a+) YE, h,,(x,0) (4.7) 
b ft 


subject to Ax =b 


x20, 
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where h,,(x,@) = minimize MOGPEN,,* Ry, 


subject to 


(1) > »y >| “oc Req,,*(GT ime_,,, + MRT ime’, ) / 24}. bee 


a reR, eT 
t'+|DTame,,, |=t 


+¥Y DY [MocReg,,*(GTime,, + MRTime’, ) /24]*Y,9 
aor Vel, 
t'+| DTime,,, |= 


— Ry < MOGCap, 
(2) Ry 20, 


where the first-stage decision vector, x, 1s used to represent all decision variables 
associated with the deterministic Throughput II model (see Appendix A). These include 
aircraft scheduling, allocation and inventory decisions as well as accounting variables to 
track deliveries of troops and equipment. 

The stochastic optimization model (4.7) 1s unique because it has (6*¢) convex 
piecewise linear cal functions. Expanding upon the theory in Section B, model 


(4.7) can be written in the following form that suggests Benders algorithm: 
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min 60 dia Y > 4,, (4.8) 
b ft 


subject to Ax =b 
cuts —G,,x+0,, > 2, VBE IA 2 clap 


where G,, is the cut gradient and g,, is the cut intercept. The calculation for the cut 
gradient and intercept expressions are contained in Appendix C. Application of Benders 
alporithna proceeds as follows (see Figure IV-1 for the algorithm): 

1. Initialization 


A constant called TOLER is initialized to a value greater than zero (typically, 
about 10°). This constant is used to compare the objective function value upper and 


lower bound, denoted Z and z respectively. Also, this step initializes @,, to zero for all 


airfields 6 and time periods ¢; and sets the objective function value upper bound 7 to a 


very large number. 
2. Master Problem 
The Master problem is solved and the optimal solution is denoted (2,6 5). since 


6, will be less than or equal to the associated recourse cost, the objective function value 


lower bound can be calculated as 
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ze + > ¥6,,. (4.9) 
5b ft 


5. Subproblem 


aA 


The first-stage decision solution x 1s passed to the subproblem. The recourse 
function is solved, obtaining the value of R,’ for all airfields 5, time periods ¢ and 
scenario @. With this information, along with the first-stage decision solution, the 


objective function value 2 can be calculated as 


Fe c&+ >) YY py MOGPEN* Ry . (4.10) 


b tf @ 


Since the upper bound = does not necessarily decrease in successive iterations, a 
comparison is made between 7 and 2; if 7<Z then the upper bound is set equal to the 


current objective function value Z and the associated solution x is saved. 
4. Optimality Check 


A comparison is made between the upper and lower bounds on the optimal 


objective function value. Specifically, if --z < TOLER*z then the algorithm 


terminates. If not, the algorithm proceeds to Step 5. 
5. Addition of Cuts 


The following set of cuts is added to the Master problem 
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-G,x +0,22 Vt (4.11) 


and the algorithm solves the Master problem (go to Step 2). Again, see Appendix C for 


the details of this calculation. 





Define TOLER = 0 






Initialize set of cuts with 6,,2>0 V br 






ZF € +00 







Solve the master problem and obtain (2,6,,] 
re cxt > »y 6, 
b ft 
Step3 Calculate R? V5, 
2 c&+¥) YY p MOGPEN* RE 


6 ft @ 











If 2 <= then =< and x* <x 






Step4 If F-z s< TOLER*z thenstop. x’ is solution. 





Augment the set of cuts with -G,x +6, 22), V bt 






Goto Step 2. 


Figure [V-1: Benders Decomposition Algorithm. 
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D. PERFORMANCE 


The algorithm was coded using the FORTRAN 77 programming language and 
IBM’s Optimization Subroutine Library (IBM, 1991). 

To improve the efficiency of the algorithm, two techniques are used. Each time 
the algorithm reaches Step 5 (see Figure IV-1), the algorithm will pass (b*1) cuts to the 


Master problem, resulting in possible redundant constraints. To remedy this potential 


problem, a cut for airfield 5 and time period f is passed to the Master problem only if 
0, < E, My (X,@) . 
The second technique is an application of Jensen’s inequality. For a fixed first- 


stage decision x the recourse functions h,,(x,@) are convex functions of the modified 


repair time vector, MRTime,,, Due to Jensen’s inequality, the expected recourse 


function £,h,.(x,@) is bounded below by the second-stage cost when we substitute 


E[MRTime,»,| for the random ground times. As stated in Chapter III, Section B.5, the 
stochastic optimization model uses the original deterministic MOG constraints. By using 
the expected value of the modified repair time random variable, in addition to the 
standard ground time, the deterministic MOG constraints provide a “more informed” 
first-stage decision in the earlier iterations of the algorithm. Then in subsequent 
iterations, cuts derived from the stochastic MOG constraints “correct” for this optimism. 
Using the algorithm shown in Figure IV-1, the stochastic optimization model 


solves in under 20 minutes on a IBM RS6000 model 590 workstation (the deterministic 


>) 





model, using GAMS, solves in under two minutes), obtaining a near-optimal solution 
(TOLER = .01). With MOGPEN,, = .04 for all airfields 6 and time periods ¢, the 
algorithm terminates in only nine iterations. As shown in Chapter III, Section C, the 
stochastic optimization model has a total of 357,591 constraints. By using Benders 
algorithm, the first iteration started with 7,671 constraints (the deterministic constraints 


of Throughput II) and the last iteration only had an additional 1,087 stochastic cuts. 
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V. SIMULATION 


This thesis describes two optimization models of the strategic airlift system, a 
deterministic and stochastic optimization model. Both models produce a recommended 
deployment schedule, based on four sets of decision variables (Xyort, Yor, Allotep, and 
Release any. See Appendix A for more details.). For example, one decision variable, X,,a;, 
will require five C5 aircraft to fly from the onload airfield Travis to the offload airfield 
Riyadh starting in period two carrying Unit Bravo. Then another decision variable, Y,,,, 
will have these five C5 aircraft fly back from Riyadh to Travis starting in period three. 
However this recommended deployment schedule was developed by an optimization 
model with its associated modeling techniques and assumptions. So, can this deployment 
schedule be executed by the Air Mobility Command of the United States Air Force? To 
address this question, a discrete-event simulation model of the strategic airlift system is 
developed. The simulation model attempts to execute a given deployment schedule and 
allows the user to analyze the recommended schedule. 

For this thesis the simulation model is used to execute the deployment schedule 
from the optimization models. The output from the simulation model is then used both 
to analyze the deterministic and stochastic models and also to compare their proposed 
deployment schedules. Just as important, examining the results of an optimization model 
via a simulation model can provide useful feedback with respect to the validity of various 


linear programming modeling assumptions. The methodology used to develop the 
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simulation model is described in Section A. Section B gives a broad overview of the 
simulation model followed by a more detailed discussion in Section C. Then in Section 


D, the performance of this simulation model is examined. 
A. METHODOLOGY 


The simulation is a discrete-event simulation model of the strategic airlift system. 
Given a deployment schedule, the simulation will attempt to execute this schedule. Two 
key concepts are employed ‘in the model. The first key concept is a “blind” execution of 
the proposed deployment schedule. For example, all aircraft will fly specific routes, even 
if that route is congested at downstream airfields and a better route is available. This 
“blind” execution will allow analysis of the proposed schedule. A properly planned 
deployment schedule will have aircraft land and immediately be serviced. On the other 
hand, a poorly planned schedule will result in delays as aircraft wait for service at 
congested airfields. 

The second key concept is the ability to execute the simulation model with 
various options. One option allows the user to turn “ON” or “OFF” the MOG constraint 
at each airfield. When the MOG constraint is “OFF,” all aircraft will land and 
immediately be serviced even when the MOG limit at that airfield is exceeded. On the 
other hand, the MOG constraint “ON” option enforces the MOG limit at each airfield. 
The ability to turn the MOG constraint “ON” or “OFF” will allow the user to examine the 
effect of a MOG constraint on the proposed schedule. A second option allows the user to 


run the simulation model with aircraft reliability “ON” or “OFF. When the aircraft 
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reliability option is “OFF,” aircraft never break upon landing; on the other hand the 
aircraft reliability “ON” option allows for the possibility for aircraft to break. This 
option will allow the user to analyze the effect of aircraft reliability on the proposed 


deployment schedule. 
B. OVERVIEW — 


This section gives a broad overview of the simulation model. A more detailed 


discussion is contained in Section C. 
1. Deployment Schedule 


The simulation model executes a given deployment schedule. Therefore the 
simulation does not need to determine which route a loaded aircraft will fly; or the 
amount of cargo to load onto an aircraft. These decisions have been determined by the 
optimization models and are represented by four types of decision variables. Therefore 
this deployment schedule given by four types of decision variables is the driving force in 


_ the simulation model. 
a. Airlift Mission Orders (X Orders) 


One type of decision variable is the actual airlift missions flown by the 
aircraft from onload to offload airfields. This decision vector [X,on] states how many 


aircraft of type a will transport unit w starting in period / using route r. To convert this 
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| variable to the simulation model, there will be_X,,,, individual X orders.? Each individual 
X order will arrive at airfield 5 in time period t. The X order requires one aircraft of type 
a to transport cargo and/or troops from unit u along route r. When the X order arrives at 
airfield 6 in time period ¢, one of two outcomes will occur. If there is an available 
aircraft of type a at airfield 6, this X order is assigned to the aircraft. Once this 
| assignment is made, the aircraft will execute the order. On the other hand if there are no 
aircraft of type a available at airfield b, the X order is assigned to a first-in, first-out 
(FIFO) queue at airfield 6. When an aircraft of type a becomes available at airfield b at a 
later time, the X order is assigned to the aircraft. 

Once the X order is assigned, the aircraft will consume MOG resources 
and onload cargo for unit w. Once the onload is complete the aircraft will take off and fly 
along route r If the required route has any enroute airfields the aircraft will land at the 
enroute airfield and the repair time random variable is used to determine the required 
repair time if the reliability option is “ON.” The aircraft will then consume MOG 
resources to refuel the aircraft and, if needed, repair the aircraft. When complete the 
aircraft will take off and ‘esdiaias to fly along route r. Eventually the aircraft will land at 
the offload airfield. Once again the repair time random variable is used to determine the 


required repair time if the reliability option is “ON.” The aircraft will consume MOG 


* The decision variables are continuous numbers. The simulation model will convert the continuous 
number to an integer by (1) rounding, (2) rounding down or (3) rounding up the solution from the 
optimization model. The simulation model allows the user to choose one of these three options for each 
type of decision variable X,,ar, Yar, Allotzn, and Release gay. 
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resources to offload the cargo and, if needed, complete repairs. When complete one of 
two outcomes will occur. If there is an old Y order (to be discussed in the next 
subsection) at airfield 5 for an aircraft of type a, the order will be assigned to the aircraft 
and the aircraft will execute the order. On the other hand if no Y order exists, the aircraft 


will remain in inventory at airfield 6 and wait for a future order to arrive. 
b. Aircraft Recovery Orders (Y Order) 


Another type of decision variable is the recovery flights of aircraft from 
the offload airfields back to the onload airfields. This decision vector [Y,,;] states how 
many aircraft of type a are to fly recovery route r starting in period ¢. To convert this 
variable to the simulation model, Y,,, individual Y orders will arrive at airfield 5 in time 
period ¢. Each individual Y order requires an aircraft of type a to fly recovery route r. 
When the Y order arrives at airfield b in time period ¢, one of two outcomes will occur. 
If there is an aircraft of type a at airfield b, this Y order is assigned to the aircraft. Once 
this assignment is made, the aircraft will execute the order. On the other hand, if there 
are no aircraft of type a available at airfield b, the Y order is assigned to a FIFO queue at 
airfield 6. When - aircraft of type a becomes available at airfield 5 at a later time, the Y 
order is assigned to the aircraft. 

Once the Y order is assigned, the aircraft will immediately take off and fly 
along route r. Eventually the aircraft will land at the onload airfield and the repair time 
random variable is used to determine the required repair time if the reliability option is 


“ON.” At this point, one of two outcomes will occur depending on whether the aircraft is 
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mission capable or not. If the aircraft is mission capable, the aircraft will either be 
assigned an X order that arrived at the airfield at an earlier time for an aircraft of type a, 
or the aircraft will be placed in inventory at airfield b and wait for an X order to arrive. 
On the other hand, if the aircraft requires repair work, the aircraft will consume MOG 
resources to complete repairs. Once the repairs are complete, the aircraft will either 
execute an X order that arrived at the airfield at an earlier time for an aircraft of type a, 


or the aircraft will be placed in inventory at airfield 5 and wait for a future X order 


to arrive. 
C, Aircraft Allotment Orders (A Order) 


The third type of decision variable is the allotment of aircraft to the airlift 
system. This decision vector [A//ot,,] states how many aircraft of type a to add to onload 
airfield 5 in time period ¢. During the execution of the simulation, A//ot,,, individual A 
orders will be executed in time period ¢. Each A order will add one aircraft of type a to 
onload airfield 6 in time period ¢. When the aircraft arrives in time period ¢ at onload 
airfield b, one of two outcomes will occur. If there is an old X order for an aircraft of 
type a at airfield 5, the order is assigned to the aircraft and the order is executed. On the 
other hand, if there are no X orders at airfield 5, the aircraft will be placed in inventory at 


airfield 6 until an X order arrives. 
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d. Aircraft Release Orders (R Order) 


The last group of decision variables releases aircraft from the airlift 
system. This decision vector [Release,,,] states how many aircraft of type a are to be 
released from onload airfield 5 in time period ¢. During the execution of the simulation, 
Release,», individual R orders will arrive at enlond airfield 5 in time period ¢. One of two 
outcomes will occur. If there exists an available aircraft of type a at airfield 5, the 
aircraft will leave the airlift system. On the other hand, if there are no aircraft of type a 
available, the order will wait in a queue (FIFO) until an aircraft of type a arrives. When 
an aircraft of type a is available, the assignment of the R order takes precedence over an 


X order. 
ym Aircraft Service Priority 


One very important modeling issue is the priority for processing aircraft at the 
airfields. When the MOG constraint option is “OFF,” aircraft will land and immediately 
consume MOG resources as required. But when the MOG constraint option is “ON,” 
some rule must be used to prioritize the aircraft. When aircraft require service (onload, 
refuel, offload, or aircraft repairs) in the simulation, the aircraft are assigned either a 
code-one or code-two service priority. The distinction is as follows. When an airfield is 
MOG constrained, code-one aircraft are always serviced first, followed by code-two 
aircraft. Also, if a code-one aircraft arrives at a MOG constrained airfield, the simulation 


will interrupt any code-two aircraft that are consuming MOG resources. For example, 
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Suppose an airfield is MOG constrained with a specified number of code-two aircraft 
consuming MOG resources. If a code-one aircraft lands, the required number of code- 
two aircraft will be interrupted to allow the code-one aircraft to consume 
MOG resources. 

The assignment of service priorities is as follows. A code-two service priority is 
always assigned to an aircraft that is non-mission capable (i.e., breaks) upon landing. 
This will allow code-one aircraft to immediately be serviced and complete their mission. 
A code-one service priority can be assigned to an aircraft in one of two ways. First, if an 
aircraft is mission capable upon landing and requires MOG resources, the aircraft is 
assigned a code-one service priority. The second way an aircraft can be assigned a code- 
one service priority is as follows. If an aircraft is non-mission capable after landing, a 
code-two service priority is assigned to the aircraft. Now if the airfield is MOG 
constrained, the aircraft will wait until all code-one aircraft have been serviced, but this 
can be an indefinite amount of time since aircraft may continue to arrive. Therefore, one 
assumption in the simulation model is that all aircraft must be repaired within the repair 
time pndonvanaBle realization. With this assumption, the simulation will upgrade a 
code-two aircraft to code-one after a predetermined amount of repair/waiting 
time elapses. 

To determine the time a code-two aircraft is upgraded to code-one is as follows. 
When an aircraft lands at an airfield and breaks, the required repair time realization and 


associated modified repair time realization are known. The repair time realization is the 
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number of hours until the aircraft is mission capable. The modified repair time 
realization is the hours of MOG resources required to repair this aircraft. As stated 
earlier, the simulation will attempt to repair all aircraft within the repair time random 
variable realization. Thus if the repair time realization is two hours and the modified 
repair time realization is two hours, this aircraft must consume MOG resources 
immediately to complete repairs. In this case the aircraft, upon breaking, is assigned a 
service priority of code-one, not code-two. On the other hand if the repair time random 
variable realization is 72 hours and the modified repair time realization is 30 hours, this 
aircraft can wait at most 42 hours prior to consuming MOG resources. In this case, the 
aircraft will transition to code-one after waiting 42 hours to consume MOG resources. 

In general the difference between the remaining repair time and the remaining 
modified repair time is called the Slack Time. When the aircraft initially breaks, the 


Slack Time is equal to the repair time realization minus the modified repair time 


realization. When the aircraft 1s waiting to consume MOG resources, the remaining 


repair time is decreasing linearly with time. As the remaining repair time decreases 
linearly with time, eventually the Slack Time will equal zero. At this time the aircraft 
will be upgraded from code-two to code-one. On the other hand if the aircraft is 
consuming MOG resources, the remaining repair time and modified repair time are both 
decreasing linearly with time. Therefore while the aircraft 1s consuming MOG resources, 


it is impossible for the aircraft to be upgraded from code-two to code-one. 
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Two criteria are used to determine when a non-mission capable aircraft becomes 
mission capable. First, the modified repair time requirement must be satisfied. This can 
only be accomplished when the aircraft consumes MOG resources for the required 
amount of time. The second criteria is the repair time must be satisfied. This is 
equivalent to requiring that the Slack Time equal zero. The second criteria is needed for 
the following reason. If a code-two aircraft immediately consumes MOG resources upon 
landing, the modified repair time criteria has been met. But the modified repair time is a 
measure of the impact of this broken aircraft on MOG. In reality, the plane will not be 


available until the repair time requirement has been satisfied. 


3. Output 


The purpose of the simulation model is to examine the feasibility of an airlift 
deployment schedule and to compare the solution provided by the deterministic and 
stochastic solutions. With this in mind, two key outputs from the simulation model are 
the delivery profile and landing profile graphs. The delivery profile graph plots the total 
amount of cargo and passengers delivered (in stons) to the offload bases versus time. 
The landing profile graph plots the total number of aircraft landings at the offload bases 
versus time. These graphs are extremely useful in examining the proposed schedule from 
an optimization model and to compare the solution from the deterministic and stochastic 
optimization models. This will be discussed further in Chapter VI. 

Another useful output from the simulation model is a graph of the number of 


aircraft that are waiting for or consuming MOG resources versus time for each airfield. 
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A quick look at the graph allows the user to determine at what time an airfield is MOG 


constrained and by how much. 

C. MODEL 
This section gives a detailed description of the simulation model. 
1. Input Data 


The required input data to run the simulation model is as follows: 


e Airfield data. To execute a deployment schedule, the only data required is the 
latitude/longitude and MOG value of each airfield. The latitude/longitude 1s 
used to compute the distance between each airfield. 


e Aircraft data. The only data required is the block speed and the deterministic 
ground time for the each aircraft type. The block speed will allow the 
simulation model to correctly compute the time for an aircraft to fly 
between airfields. 

e Aircraft reliability. This input file gives the repair time random variable 
realizations, modified repair times realizations and associated probabilities for 


C5, C17 and C141 aircraft. All other aircraft do not break. 


e Routes. This data is used in conjunction with the optimization solution so the 
simulation can properly route aircraft from airfield to airfield. 


e Optimization solution. The driving force for the simulation is the deployment 
schedule developed by the optimization models. This consist of four groups 
of decision variables Xjart, Yar, Allotgn,, Release gp. | 


Z; Queuing System for MOG Consumption 


There are a total of eight queues used to properly process aircraft. They are 


described below and a diagram 1s shown in Figure V-1. 
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a. Taxi Queue 


If the airfield is MOG constrained, aircraft are placed in the Taxi (FIFO) 
Queue. When a unit or partial unit of MOG becomes available, the next aircraft in the 
Taxi Queue will enter the Code-1F, Code-1P, Code-2F or Code-2P Queue. If the next 
aircraft in line has a service priority of code-two, the aircraft immediately enters the 


Hold-In Queue and time now counts toward the repair time requirement. 
b.  Code-1F Queue 


This queue is for aircraft with a service priority of code-one. Each aircraft 
is consuming a full unit of MOG and time is counting toward both the modified repair 
time and repair time requirement (if any). The aircraft will remain in this queue for an 
amount of time equal to the modified repair time plus the standard ground time. When 


services are complete, the aircraft will then exit. 
c. Code-1P Queue 


This queue is for aircraft with a service priority of code-one, but the 
aircraft is only consuming a partial unit of MOG. At any given time, there will be at 
most one aircraft in this queue. Time is counting toward both the modified repair time 
and repair time requirement (if any). If an aircraft leaves the Code-1F Queue at a MOG 
constrained airfield, the aircraft in the Code-1P Queue, if it exists, will enter the Code-1F 


Queue and consume a full unit of MOG. It is also possible for an aircraft to remain in the 
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Code-1P Queue long enough to satisfy the modified repair time requirement plus the 


standard ground time, at which time it will exit. 
d. Code-2F Queue 


This queue is for aircraft with a service priority of miewe where each 
aircraft is consuming a full unit of MOG. Time is counting toward both the modified 
repair time and repair time requirement. The aircraft will remain in this queue for an 
amount of time equal to the modified repair time plus the standard ground time. When 
services are complete the aircraft will enter the Hold-Out Queue until the repair time 
requirement is satisfied. If at any time a code-one aircraft needs to consume resources 
and the airfield is MOG constrained, the required number of aircraft in the Code-2F 


Queue, if any, will be interrupted and enter the Hold-In Queue. 
é. Code-2P Queue 


This queue is for an aircraft with a service priority of code-two but this 
aircraft is only consuming a partial unit of MOG. At any given time, there will be at 
most one aircraft in this queue. Time is counting toward both the modified repair time 
and repair time requirement. When an aircraft leaves the Code-2F or Code-1F Queue and 
the airfield is MOG constrained, the aircraft in the Code-2P Queue, if it exists, will enter 
the Code-2F Queue and consume a full unit of MOG. It is also possible for the aircraft to 
satisfy the modified repair time requirement plus the standard ground time while the 


aircraft is in this queue, at which time the aircraft will enter the Hold-Out Queue until the 
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repair time requirement is satisfied. If at any time a code-one aircraft needs to consume 
resources and the airfield is MOG constrained, the aircraft in the Code-2P, if there is one, 


will be interrupted and enter the Hold-In Queue. 
f Hold-In Queue 


Only code-two aircraft enter this queue. Time is counting toward the 
repair time requirement. When a unit or partial unit of MOG becomes available in either 
the Code-2F or Code-2P Queue, an aircraft, if any, will leave the Hold-In Queue and 
enter the appropriate queue. If the Slack Time reaches zero for any aircraft while in this 


queue, the aircraft will be upgraded to code-one and enter the Upgrade Queue. 


2. Upgrade Queue 


Only aircraft that have been upgraded from code-two to code-one enter 
this queue. From the Upgrade Queue aircraft will enter either the Code-1F or Code-1P 


Queue when a unit or partial unit of MOG becomes available. 
h. Hold-Out Queue 


This queue is for code-two aircraft that have satisfied the modified repair 
time requirement, but have yet to satisfy the repair time requirement. The aircraft is not 
consuming MOG resources. When the repair time requirement has been satisfied, the 


aircraft will exit. 
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Hold-Out Queue 


Code-1F Queue _ 
{consume MOG) | 


Code-1P Queue | 
~ (consume MOG) — 


Code-2F Queue 
—*_tersane os 


Code-2P Queue 
(consume MOG) 


Upgrade Queue 
Hold-In Queue 


Time counts toward the repair time requirement. 





Figure V-1: Queuing System for MOG Consumption. 
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i. 





Priorities Among Queues 


When an aircraft has completed consuming MOG resources, the 


simulation will use the following priority to determine the next aircraft to consume 


MOG resources: 


ae 2S he 


Code-1P Queue 
Code-2P Queue 
Upgrade Queue 
Taxi Queue 

Hold-In Queue. 


3. State Variables 


In simulation models, state variables are a collection of variables necessary to 


describe a system at a particular time. For this simulation model the state variables can 


be grouped into two general categories. 


a. 
Aircraft Type 
ID Number 


Cargo 


Cargo ID 


Break Code 


Status 


Aircraft State Variables 


Defines the aircraft type; e.g., C5, C17, C130, C141, 747C, 
747P, DC10. 


The identification number of the aircraft. This allows the 
simulation to track each individual aircraft. 


This variable is true if the aircraft has cargo onboard and fully 
loaded, and false otherwise. 


The name of the unit being transported by the aircraft. 


Defines the required repair time for the aircraft (MC, two, six, 
ten, fourteen, twenty, thirty-six, sixty, and seventy-two). 


Defines the stage of the aircraft in the system. They are as 


follow: airborne, onload, refuel, offload, repair, inventory, pre- 
airlift (the aircraft has not yet entered the airlift system) and 
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post-airlift (the aircraft has been released from the 
airlift system). 


Service Priority When an aircraft needs to consume MOG resources, this state 
variable gives the location of the aircraft in the queuing system 
for MOG consumption (see Figure V-1). The possible values are 
taxi, code-1F, code-1P, code-2F, code-2P, hold-in, hold-out, 
upgrade and NA (the aircraft does not require MOG services). 


Current Base The name of the airfield where the aircraft is currently located. 


Entered Service The time the aircraft left the Taxi Queue and entered the Code- 
1F, Code-1P, Code-2F, Code-2P or the Hold-In Queue. 


MOG Time The total hours of MOG resources required by this aircraft. This 
includes the standard ground time and the modified repair time. 


Cum Time The current amount of MOG resources consumed by the aircraft. 


b. Airfield State Variables 


Name The name of the airfield. 

Lat The latitude position of the airfield. 

Long The longitude position of the airfield. 

MOGReq The narrow-body equivalent requirement for each aircraft type 
for this airfield. 

MOG The agtrow-body MOG limit. 

Current The number of narrow-body equivalent aircraft consuming 


MOG resources. 


Req'd The number of narrow-body equivalent aircraft that are 
consuming MOG resources and those aircraft that need to 
consume MOG resources. Specifically, this is the narrow body 
equivalent sum of aircraft whose service priority state variable is 
equal to code-1F, code-1P, code-2F, code-2P, taxi, and upgrade. 
This is used as a measure of how much the MOG value is 
being exceeded. 
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R Orders _ There is a FIFO queue for each aircraft type. Each queue holds 
old R orders that still need to be executed. 

X Order There is a FIFO queue for each aircraft type. Each queue holds 
old X orders that still need to be executed. 

Y Order There is a FIFO queue for each aircraft type. Each queue higldé 
old Y orders that still need to be executed. 

Inv Lists all aircraft that are being held as inventory at this airfield. 

Taxi Lists the aircraft assigned to the Taxi Queue for that airfield. 

CodelF Lists the aircraft assigned to the Code-1F Queue for that airfield. 

Code1P Lists the aircraft assigned to the Code-1P Queue for that airfield. 

Code2F Lists the aircraft assigned to the Code-2F Queue for that airfield. 

Code2P Lists the aircraft assigned to the Code-2P Queue for that airfield. 

Upgrade Lists the aircraft assigned to the Upgrade Queue for that airfield. 

Holdin Lists the aircraft assigned to the Hold-In Queue for that airfield. 

HoldOut Lists the aircraft assigned to the NHold-Out Queue for 
that airfield. 

4. Events 


In simulation modeling events are defined as an instantaneous occurrence that 
may change the state of the system. The events used in this simulation model are briefly 
described below and the event diagrams are contained in Appendix E. 


This event adds an aircraft to an onload airfield. The aircraft 
will either remain in inventory at the airfield or be assigned an 
old X order. If an old X order exists the aircraft will execute that 
X order. The event diagram 1s contained tn Figure D-1. 


Allot Order 


70 





Release Order 


X Order 


Y Order 


Land 


Ready-1F 


~Ready-1P 


Ready-2F 


Ready-2P 





This event attempts to remove an aircraft from the airlift system. 
If the aircraft is not available, the order will enter a FIFO queue 
at the airfield until the appropriate aircraft arrives. The event 
diagram is contained in Figure D-2. 


This event attempts to execute a new X order. If there is an 
available aircraft, the X order is assigned to the aircraft and the 
aircraft will execute the X order. If there are no aircraft 
available, the X order is entered in a FIFO queue at the specified 
airfield. The event diagram is contained in Figure D-3. 


This event attempts to execute a new Y order. If there is an 
available aircraft, the Y order is assigned to the aircraft and the 
aircraft immediately takes off. On the other hand, if there are no 
available aircraft, the Y order is entered in a FIFO queue at the 
specified airfield. The event diagram is contained in Figure D-4. 


This event lands an aircraft at an airfield. If the reliability option 
is “ON,” the repair time random variable is used to determine the 
amount of time required to repair the aircraft. If the aircraft just 
arrived at an onload airfield and is mission capable, the aircraft 
is available to execute an X order. If the aircraft needs to refuel, 
offload or complete repairs, the aircraft will consume MOG 
resources. The event diagram is contained in Figure D-5. 


A specific aircraft in a Code-1F Queue just completed 
consuming MOG resources and will exit. The aircraft will 
continue to execute its assigned order. The event diagram is 
contained in Figure D-7. 


A specific aircraft in a Code-1P Queue just completed 
consuming MOG resources and will exit. The aircraft will 
continue to execute its assigned order. The event diagram is 
contained in Figure D-7. 


A specific aircraft in a Code-2F Queue just completed 
consuming MOG resources. The aircraft will enter the Hold-Out 
Queue. The event diagram is contained in Figure D-8. 


A specific aircraft in a Code-2P Queue just completed 
consuming MOG resources. The aircraft will enter the Hold-Out 
Queue. The event diagram is contained in Figure D-8. 
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Ready-Hold —_A specific aircraft in the Hold-Out Queue has satisfied the repair 
time requirement and will continue its mission. The event 
diagram is contained is Figure D-10. 

Transition A specific aircraft in the Hold-In Queue is upgraded from a 
service priority of code-two to code-one. The aircraft is removed 


from the Hold-In Queue and entered in the Upgrade Queue. The 
event diagram is contained in Figure D-11. | 


5. Execution Options 


There are numerous options available in the simulation. The main options are 


described below. 
a. Time to Execute Orders 


As stated earlier, the optimization models use discrete time periods, but 
the simulation model has a continuous time spectrum. For example, an optimization 
model solution might add four CS aircraft to airfield 5 in day two. But in the simulation, 
day two extends from 24 to 48 aiaes At what time should the simulation model add 
these aircraft? To address this issue the simulation model provides the user with two 
options. Using the first option, the simulation model will have all new orders arrive at 
the beginning of the corresponding time period. Using the second option, the simulation 
model will uniformly spread the new orders throughout the time period. For adialie 
suppose the optimization model requires that four C5 aircraft be added to airfield 4 in 


day two. Using the second option, the simulation will add the four C5 aircraft at hours 


24, 30, 36 and 42. 
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b. Converting the Continuous Decision Variable Solutions to 


Integers 


The four groups of decision variables are used as input data for the 
simulation model. But the variables are continuous numbers and the simulation must use 
integer numbers. Therefore for each one of the four decision variables, the user has three 
options. They are (1) round-down, (2) round, and (3) round-up. This will be discussed | 


further in Chapter VI. 
Cc. MOG constraint 


If this option is “ON” the simulation model will enforce the MOG 
constraint at each airfield. On the other hand if this option is “OFF,” all aircraft will 
consume MOG resources as needed even when the MOG limit at that airfield 


is exceeded. 
d. Aircraft Reliability 


When the aircraft reliability option is “OFF,” aircraft never break upon 
landing; on the other hand the aircraft reliability “ON” option allows for the possibility 


for aircraft to break. 
D. PERFORMANCE 


The simulation model is coded using Borland Pascal with Objects (Borland, 
1992). If the simulation model is run with the reliability option “OFF,” only one run is 


required and this run takes about 12 seconds (486 DX2 66 Hz computer). If the 
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simulation model is run with the reliability option “ON,” multiple runs are required until 


Steady state is achieved and this takes about six hours. 
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VI. ANALYSIS 


This chapter reports the results of using the simulation model as a tool to saves 
the deterministic and stochastic optimization models. For all analyses a base scenario 
Meveiepea by the U.S. Air Force Studies and Analyses Agency is used. Section A 
discusses two separate issues that must be resolved before the simulation can be used as 
an analytical tool. Then in Section B the simulation model is used to examine the effect 
of the deterministic MOG efficiency factor. Section C reports the results of using the 
simulation model to analyze the deterministic model. Finally, Section D compares the 


deterministic and stochastic optimization models using the simulation. 

A. COUPLING OF THE OPTIMIZATION MODELS WITH SIMULATION 
Scctiog A explains default settings for two options that may be set by the user. 
1. Converting the Continuous Decision Variable Solution to Integers 


As stated earlier, the four types of decision variables passed to the simulation 
from the optimization model are continuous numbers. To convert the continuous 
decision variable solution to integers, the simulation model wil! round-down, round, or 
round-up depending on the option chosen by the user. The following criteria is used to 
choose the correct setting for this option for each type of decision variable. Since the 
simulation model is used to study the deployment schedule recommended by the 


optimization models, the simulation model needs to execute the same number of airlift 
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missions (X orders) and recovery flights (Y orders). For example if the total number of 
airlift missions in the deterministic model is 100.40, the simulation model needs to 
execute a total of 100 or 101 X orders. Once this is achieved, there will still be an error 
associated with the individual rounding of the decision variable solution. For example, 
one airlift mission might require 12.5 C5 aircraft to airlift unit w using route 7 starting in 
time period ¢, but if the round-up option is used, this will be translated to 13 C5 aircraft. 
To examine the impact of this issue, the deterministic optimization model is 
solved with a MOG efficiency factor equal to one. Using this solution as input, the 
simulation model is executed using the three different round options for each one of the 
four types of decision variables. When the round-down option is used for the X orders, a 
total of 3,624 X orders are generated in the simulation model compared to the original 
number of 3,734.40 from the optimization model. In this case the simulation will 
attempt to execute 3,624 X orders and is 110.40 short of the required amount. Using the 
round-down option for the Y orders, the simulation attempts to execute 3,280 Y orders 
compared to 3,388.70 from the optimization model. In this case the Y orders are short by 
108.7. In summary ieeqauna dean option results in a total of 219.1 X and Y orders that 
are not considered for execution in the simulation model. The same analysis can be 
performed for the round and round-up option. See Table VI-1 for a summary of 


the results. 
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a Bs Round-Down ; | Round-Up 
Sum of X Orders - Simulation 3/34 
Sum of X Orders - Optimization 3734.4 
X Order Excess ; -0.4 


Sum of Y Orders - Simulation 3387 3499 
Sum of Y Orders - Optimization ; 3388.7 3388.7 
Y Order Excess -1.7 110.3 


oe 
219.1 
aa: 
12 
‘ 


Number of X Orders Not Executed 
Number of Y Orders Not Executed 
Number of X and Y Orders Not Executed 





Corrected X and Y Order Excess -241.1 


Table VI-1: Simulation Performance versus Round Option. The deterministic 
optimization model (MOGEff = 1.00) Is used as input into the 
simulation model. 


There is a consequence of converting the continuous decision variable solution to 
integers. The optimization models use aircraft balance constraints at onload and offload 
airfields to ensure there are aircraft available to fly the airlift missions and recovery 
flights. When the decision variable solution is converted from a continuous number to an 
integer, this aircraft balance might no longer exist. In the simulation, for example, an X 
order will not be executed if there was no corresponding Y order to return the required 
aircraft type to that airfield. Therefore at the completion of one run, the simulation will 
count the number of X and Y orders that are not executed. The results are also 
summarized in Table VI-1. 

As shown in Table VI-1, when the round-down option is used for X and Y orders 


the simulation model needs to execute an additional 241.10 X and Y orders to properly 


ed 





model the flow of aircraft through the airlift infrastructure. The round-up option 
executes an extra 192.90 X and Y orders. However, by using the round option, the 
simulation is only short by 7.10 orders. Therefore when using a MOG efficiency factor 
equal to one, the round option closely matches the optimization model. Similar results 
are obtained if different MOG efficiency factors are used. Therefore all analyses use the 
round option for X and Y orders. 

In converting the continuous decision variable solution associated with the 
allotment and release of aircraft, all analyses use the round option. In most cases there is 
no difference between the number of aircraft used in the optimization and the 


simulation models. 
=? Arrival of Orders 


A second simulation option that the user may control regards scheduling the 
arrival of orders within the time periods given by the optimization models. As stated 
earlier the optimization models use discrete time periods, but the simulation model has a 
continuous time spectrum. To address this issue, the simulation model gives the user two 
options. The first option will have all oe arrive at the very beginning of the time 
period, denoted “0000.” The second option will have the orders arrive uniformly 
throughout the trme period, denoted “uniform.” Table VI-2 illustrates the option used in 


all analyses for each decision vartable. 
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Release apt uniform 


Table VI-2: Arrival of Orders Option. 







The reason for this choice is as follows. The simulation should use a realistic 
schedule within the time periods. In airlift operations the Air Force attempts to maintain 
sufficient spread between departure times so as to minimize congestion at onload, 
enroute and offload airfields. This same type of concept is applied in the simulation 
model by using the “uniform” option for A (If there are old X orders at the airfield, 
aircraft will be assigned the X orders and request MOG services to onload the cargo.) and 
X orders. On the other hand, when the aircraft arrive at the destination and offload their 
cargo in airlift operations, they will immediately take off and return to an onload airfield. 
Therefore the “0000” option is used for the arrival of Y orders at offload airfields. 
Finally, the “uniform” option is used for the arrival of R orders to uniformly spread the 
departure of aircraft from the airlift system throughput the time period. 

The opnons shown in Table VI-2 result in the simulation behaving in a consistent 
fashion. Using the simulation model with the MOG constraint option “OFF” and the 
reliability option “OFF,” all aircraft will immediately consume MOG as required and 
experience no delays caused by breakage. This results in a verbatim execution of the 
airlift schedule. By examining simulation output graphs similar to Figure VI-1 for all 


airfields, analysts can determine when the given schedule is exceeding the capacity of an 
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airfield. By varying the arrival options, the user can determine the best combination for 
all decision variables. For example the “0000” option for all four decision variables 
yields the solid line in Figure VI-1. In this case the capacity of Travis airfield is 
exceeded on days six, eight and eleven. Upon closer examination of the scenario, 
aircraft are allotted to Travis airfield on day six and eleven at midnight. But there are old 
X orders existing at the airfield. Therefore the X orders are immediately assigned to the 
aircraft and the capacity of the airfield is exceeded as the aircraft onload their cargo. 
Thus the solid line in Figure VI-1 will result in delays as aircraft wait for service when 
the MOG constraint is enforced. On the other hand, the dash line is obtained by using the 
recommended options and results in a more uniform arrival rate and consumption 


of MOG. 
B. DETERMINISTIC MOG EFFICIENCY FACTOR 


As stated in Chapter II, aircraft turnaround times for onloading and offloading 
cargo and enroute refueling are assumed to be known constants in the Throughput II 
model, although they are naturally stochastic. This ignores the fact that deviations from 
the given service time can cause congestion on the ground. To offset the optimism of 
this assumption, Throughput IJ uses an efficiency factor (<1) in the formulation of MOG 
consumption constraints to lessen the impact of randomness. This section estimates a 


reasonable value for this MOG efficiency factor. 
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Figure VI-1: Simulation Output Graph of Travis Airfield. The deterministic model 
(MOGEff = 1.00) is used as input into the simulation model. The 
simulation is then run with the following options: MOG “OFF,” reliability 
“OFF.” 


1 Methodology 


The following approach is used to estimate a reasonable value for the deterministic MOG 


efficiency factor. 
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a. Optimization Model 


Solve the deterministic optimization model with a specified MOG 


efficiency factor and obtain the delivery profile graph (Figure VI-2). 
b. Effect of Discrete Time Periods 


Using the deterministic solution as input, run the simulation model with 
the following options: reliability “OFF” and MOG “OFF.” The delivery profile graph | 
obtained from this simulation run will differ from step (a); the difference between the 
two delivery profile graphs shows the effect of the discrete time periods of the 
optimization model. Examining Figure VI-2, the simulation completed the delivery of 


cargo on day 3] compared to day 26 for the optimization model. 
Cc. Effect of Aircraft Reliability 


Perform another simulation run with the following options: reliability 
“ON” and MOG “OFF.” This delivery profile graph takes into account aircraft reliability 
and reflects the minimum delay (since the MOG constraint is “OFF” to deliver all the 


cargo as aircraft spend more time on the ground to complete repairs (see Figure VI-2). 
d. Effect of MOG Constraint 


Perform a simulation run with the following options: reliability “ON” and 
MOG “ON.” The delivery profile graph now takes into account the uneven arrival of 
aircraft and the limited resources available to repair the non-mission capable aircraft. 


This effect can be adjusted by changing the MOG efficiency factor. As shown is Figure 


82 








VI-2, the end result is that all cargo will not be delivered until about day 43. One 
interesting observation is that on day 26, the optimization model has delivered all the 


cargo, but the simulation still has 35,000 stons to deliver, or approximately 17 percent. 


MOG Eff = 1.00 
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Figure VI-2: Simulation Delivery Profile Comparison. The deterministic model 
(MOGEff = 1.00) is used as input into the simulation model. 
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Li Results 


Using this approach, Appendix F contains six similar graphs with MOG efficiency 
factors of 0.50, 0.60, 0.70, .080, 0.90 and 1.00. The effect of the MOG constraint is 
negligible when a MOG efficiency factor of 0.50 is used. But operating the airfields at 
half capacity 1s unrealistic because the airlift system would be dramatically underutilized. 
On the other hand, the effect of the MOG constraint is significant when a MOG 
efficiency factor of 1.00 is used. By examining the graphs in Appendix F, one can 
conclude that the best MOG efficiency factor is approximately 0.70 or 0.80. Also if the 
six delivery profile graphs (MOG “ON” and Reliability “ON’’) are compared, the graphs 
corresponding to MOG efficiency factors of 0.70 and 0.80 deliver more cargo when the 
airlift infrastructure is MOG constrained. 

Another output from the simulation model is the landing profile at offload 
airfields versus time. This graph slots the total number of landings at the offload 
airfields versus time. Using the same approach as the delivery profile graphs, the landing 
profile graphs in Appendix D are obtained. Once again, the best MOG efficiency factor 


to use is either 0.70 or 0.80. 
C. ANALYSIS OF DETERMINISTIC OPTIMIZATION MODEL 


The assumptions and limitations of the deterministic optimization model were 
discussed in Chapter II, Section B. One limitation of the model is the discrete time 


periods. Specifically, the model enforces the MOG constraint for each airfield using the 


84 





discrete time periods. With one-day time periods, the model can route aircraft in a 
manner that causes local congestion. For example, all aircraft could arrive at an airfield 
within a small time window instead of being dispersed over an entire day. In reality this 
would cause local congestion, even though the model’s representation of aircraft 
handling capacity is observed. Another limitation is that the model rounds the time (to 
the nearest day) at which an aircraft arrives at origin and destination airfields. 

To examine this issue of discrete time periods, the simulation model can be used. 
One of the simulation output graphs is the number of aircraft that are consuming MOG 
resources and those aircraft that need to consume MOG (specifically those aircraft in the 
Taxi and Upgrade Queue) resources versus time for each airfield. Two steps are used. 
First perform a simulation run with the reliability option “OFF” and the MOG constraint 
option “OFF” and obtain the output graph. Examining Figure VI-3 the following points 
are illustrated: 


e Uneven aircraft arrival. The graph indicates a very non-uniform aircraft 
arrival rate at Mildenhall airfield. 


e MOG constraint. During the time period from day 10 to 25, the airfield is 
MOG constrained according to the deterministic optimization model solution 
report. But from day 11! to 12 the simulation model exceeds the one day 
MOG limit, but the time period from day 12 to 13 does not. This reflects the 
error caused by rounding in the optimization model. 


Next perform a simulation run with the reliability option “OFF” and the MOG constraint 


option “ON.” Now the simulation model will continuously enforce the MOG constraint, 
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unlike the one-day average enforcement by the optimization model. Examining Figure 


VI-3 the following points are illustrated: 


Up to day seven, both simulation runs track. But then they diverge even 
though the airfield is not MOG constrained. This reflects the fact that other 
airfields are MOG constrained, thus delaying the arrival of aircraft to 
Mildenhall airfield. 


From day 23 to 26 the airfield is exceeding the MOG capability of the airfield. 
At one point about 80 NB eens aircraft are on the ground waiting and 
consuming MOG resources. 


Another simulation output that is used to examine the deterministic optimization 


model is the delivery profile graph. The following steps are taken (see Figure VI-4): 


I. 


Obtain the delivery profile graph from the deterministic optimization model. 


2. Using the deterministic model solution as input into the simulation, perform a 


simulation run with the following options: reliability “OFF” and MOG “OFF.” 
This difference in the delivery profile graphs illustrates the effect of the 
discrete time periods. In this case, the graph indicates that, on average, there 
is a difference of about 10,000 stons. 


Perform a simulation run with the following options: reliability “OFF” and 
MOG “ON.” The difference in the delivery profile graphs illustrate the effect 
of limited MOG resources at the various airfields. In this case, the difference 
is negligible. 


Perform a simulation run with the following options: reliability “ON” and 
MOG “ON.” The difference in the delivery profile graph illustrates two 
effects: the delay caused by aircraft actually breaking and spending more time 
on the ground and the delay, if any, as the aircraft wait for the extra MOG 
resources required to fix the aircraft. This additional delivery profile graph 
can not be used to justify the proper choice of the MOG efficiency factor 
since the delay caused by aircraft simply breaking is not related to the MOG 
efficiency factor. 
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Figure V]I-3: Simulation Output Graph of Mildenhall Airfield. The deterministic model 
(MOGLff = 0.80) is used as input into the simulation model. The solid line 
is the simulation run with REL “OFF” and MOG “OFF.” The dash line ts 
the simulation run with REL “OFF” and MOG “ON.” 


87 








2000 


1800 


1600 


1400 


1200 


1000 


800 


600 


Passengers and Cargo (100 stons) 


400 





1 3 5 7 9 11 13 15 17 19 21 23 25 27 29 31 33 35 37 
Days 


DETERMINISTIC MODEL 
— —-— SIMULATION - REL ‘OFF,’ MOG 'OFF'} 


SIMULATION - REL ‘OFF,’ MOG 'ON' 
——— SIMULATION - REL 'ON,’ MOG 'ON' 





Figure VI-4: Simulation Deliver Profile Comparison. The deterministic model (VOGEff 
= 0.80) 1s used as input into the simulation model. 
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D. DETERMINISTIC AND STOCHASTIC OPTIMIZATION COMPARISON 


Before comparing the two optimization models, the impact of the MOG penalty 
in the stochastic optimization model is examined. Using the formulation from Chapter 


Ill, Section D, the MOG penalty can range from 0.001 to 0.26. The total amount of cargo 


delivered as a function of the MOG penalty is illustrated in Figure VI-5. 
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Figure VI-5: Total Cargo and Passengers Delivered by the Stochastic Optimization 
Model. The deterministic model is solved with a MOG efficiency factor 
equal to one. 
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To compare the stochastic and deterministic optimization models the following 
approach is used. Specify a value for the stochastic MOG penalty and solve the 
stochastic optimization model. The solution from the stochastic optimization model, 
specifically the deployment schedule, 1s used as input into the simulation model. The 
simulation model is then run with the following options: MOG “ON” and reliability 
“ON.” When the simulation has reached steady state the delivery profile graph cl 
obtained. Using this approach the three delivery profile graphs are obtained using MOG 
penalties of 0.00, 0.02 and 0.04 for the stochastic optimization model (see Figure VI-6). 
For comparison the same approach ts used for obtaining the delivery profile graph for the 
deterministic optimization model. 

The stochastic optimization model with a MOG penalty equal to zero is 
equivalent to solving a modified version of the deterministic model. The modified 
deterministic model uses the expected value of the modified repair time, in addition to - 
the deterministic ground time, for the calculation of MOG consumption. Since the 
modified repair time 1s being used tn the calculation of the MOG consumption, the MOG 


efficiency factor needs to be set to one. 
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Figure VI-6: Simulation Delivery Profile Comparison. The four delivery profile graphs 
are obtained by first solving the corresponding optimization model. Then 
the solution from the optimization model is used as input into the _ 
simulation model. The simulation is then run with the following options: 
MOG “ON” and reliability “ON.” 


All four optimization/simulation model runs can be characterized as having three 
distinct phases as observed in Figure VI-6. In the first 12 days of the base scenario, the 


system is airframe constrained and all four model runs are very similar. But from day 12 
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to about day 20, the airlift infrastructure is the limiting factor, specifically airfield 
capacity. During this fate period all three stochastic models exceed the deterministic 
model in the delivery of cargo and passengers. From day 21 and on, the system 1s in a 
sustainment phase with diminished demands. Neither airframes nor airfield capacities 
are critical resources. One interesting point illustrated by the graphs is that the schedules 
generated by all three stochastic models delivered all the cargo in the simulation model 
by day 30; but the eine model took an 10 extra days, or an additional pneahird 
of the required 30 day time span. 

Closer examination of the time period from day 12 to day 20 lends some insight 
into the power of the stochastic optimization model. When the deterministic 
optimization model is solved, the model is operating critical airfields at maximum 
capacity during this time period (day 12 to 20). This will allow the optimization model 
to deliver the largest amount of cargo. Then the solution from the deterministic 
optimization model is used as input into the simulation and the simulation model is run 
with the MOG constraint option “ON” and the reliability option “ON.” But when the 
simulation gets to this time period all it takes is for one aircraft to break at a MOG 
constrained airfield to have a dramatic impact. This aircraft will need to consume a 
larger amount of MOG resources resulting in follow-on aircraft waiting on the ground 
while this aircraft 1s being repaired. As more and more aircraft break, the problem just 


multiplies as more and more aircraft wait on the ground for MOG resources. 
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On the other hand the stochastic optimization model takes a different approach as 
illustrated by the following three points. First, when the stochastic optimization model is 
solved, the model will not operate these critical airfields at maximum capacity because 
the associated recourse costs outweigh the benefits. Therefore the stochastic model will 
operate the critical airfields at some percentage below their maximum capacity. This 
will allow extra MOG resources to be available to repair aircraft as they break. Second, 
since the stochastic optimization model will be delivering less cargo during this time 
period it will be more selective in determining which units to move. Finally the 
stochastic optimization model is re-routing aircraft to minimize the number of unreliable 
aircraft that are using capacity limited airfields. Therefore when the sails 
optimization model solution is used as input into the simulation, the simulation does not 
experience major bottle-necks at these critical airfields and will exceed the deterministic 
model in the amount of cargo delivered. 

As seen in Figure VI-6 the benefits of the stochastic optimization model take 
effect when the infrastructure is MOG constrained. In this case, this is from day 12 to 
about day 20 (see Figure VI-7). During this time period, all three stochastic model runs 
deliver more cargo than the deterministic model (no matter what number is used for the 
MOG efficiency factor, the deterministic model deliver less cargo than the stochastic 
model). Figure VI-7 indicates that as the MOG penalty gets larger the stochastic model 
delivers less cargo during this time period. But consider the following points: 

e As the MOG penalty is increased from 0.00, the delivery profile initially drops 


off, but by day 17 all three graphs are nearly the same. At this point the 
diminishing demand begins to take effect and all three graphs approach their 
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final delivery amount. If the airlift infrastructure was MOG constrained over 
a longer time period, it is possible that a difference might be noted. 


e Using the delivery profile graph as a measure of effectiveness does not 
directly reflect the value of the stochastic optimization model. Another 
possible measure would be the probability that an airfield exceeds its capacity 
in a given time period. 
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Figure VI-7: Simulation Delivery Profile Graph. The four delivery profile graphs are 
obtained by first solving the corresponding optimization model. Then the 
solution from the optimization model is used as input into the simulation 
model. The simulation is then run with the following options: MOG “ON” 
and reliability “ON.” 
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Another difference is found between the stochastic and deterministic optimization 


solutions. Examining the time period from day 12 to day 20, Table VI-3 shows the 


percent of X orders (airlift missions) flown by each aircraft type. 


Model | Ct C141 
| 
stochastic (pen = 0.00) 


stochastic (pen = 0.04) 26.3 8./ 10.7 38.4 8.5 7.4 















P747 | C/47 


C130 









Table VI-3: Percent of X orders (airlift missions) flown by each aircraft type 


from day 12 to day 20. 


The following points explain the results obtained in Table VI-3. 


Break rate data is not available for the C130, P747, and C747 and therefore were 
modeled as perfectly reliable aircraft. Relative to this, we have considerable 
insight as to why the models behave the way they do. 


The stochastic model (MOG penalty = 0.00) differs from the deterministic model 
only in that the ground times of the C5, C17, and C141 are increased by the 
expected value of the modified repair time random variable. Therefore the C130, 
P747 and C747 are, relatively speaking, more attractive in the stochastic model 
(MOG penalty = 0.00); this is why their usage increases. From Jensen's inequality 
the expected value will provide an optimistic prediction of an aircraft's 
performance relative to using the probability distribution. Thus the C130, P747, 
and C747 are again, relatively speaking, more attractive when we move from the 
stochastic model (MOG penalty = 0) to the stochastic model (MOG penalty = 
0.04); this is why the usage of the perfectly reliable aircraft again increases. 


Why does the C130 appear to displace the C141 as opposed to the C17 or C5 (the 
other two aircraft with stochastic ground times)? The C130 and C141 can both 
carry bulk and over but not out-sized cargo. To move the outsized cargo we need 
the C5 and C17. 


The C130 is a very slow (block speed) aircraft and is not usually even regarded as 
a strategic airlifter. The results do not suggest that the C130 is the answer to the 
strategic airlift problem. There are a couple of unique characteristics of this base 
scenario. First, there is the perfectly reliable assumption for the C130 aircraft. 
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Second, a large fraction of the total movement in this 30 day scenario is from 
Ramstein to Saudi Arabia; the fact that this is (in strategic airlift terms) such a 
Short distance makes the C130 more of a player than it might otherwise be. 
Finally, the fact that the theater is operating at MOG capacity means that the slow 
block speed does not slow down throughput as much as it might otherwise. 


In summary the stochastic optimization model, in addition to the features of 
Throughput I, accomplishes the following: | 
e Selection of aircraft routes by anticipating potential bottlenecks in the system. 


e Minimizes the number of unreliable aircraft that use capacity limited airfields 
or airfields that have limited repair capability. 


e Operating critical airfields at a percentage below their maximum capacity in 
order to have excess MOG resources available to repair aircraft as they break. 
This allows the stochastic optimization model to achieve a flow of cargo to 
the theater that is not interrupted by the random events of aircraft reliability. 
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Vil. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 


A stochastic extension has successfully been added to the Throughput II model to 
account for aircraft reliability. This stochastic extension results in a two-stage stochastic 
optimization model with recourse. 

The stochastic optimization model can be used to gain broad insights into the 
strategic airlift system. Some of the broad insights that can be gained by including 
aircraft reliability are the change in fleet mix and the re-routing of aircraft compared to 
the deterministic optimization model. But most important of all 1s the fast turn-around of 
the stochastic optimization model. This allows quick answers without tying down a 
massive amount of manpower or computer resources. 

In addition the airlift infrastructure has also been modeled by a discrete-event 
simulation model. This allows the user to analyze the recommended deployment 
schedule from the deterministic and stochastic optimization models. The simulation will 
also allow the user to analyze a given deployment schedule and examine the effects of 


MOG and aircraft reliability. 
B. RECOMMENDATIONS 


The following are some recommendations based on the analysis completed in 


Chapter VI. 
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1. MOG Efficiency Factor 


Using the simulation model to analyze the deterministic model, Chapter VI 


estimated a reasonable value of 0.80 for the MOG efficiency factor. 
2. Modification to Deterministic Optimization Model 


The current form of the deterministic optimization model is overly optimistic. 
Even though the MOG efficiency factor can be used to lessen the impact of random 
ground times, it does not take into account aircraft reliability. To account for aircraft 
reliability, the deterministic model can be modified by using the expected value of the 
modified repair time, in addition to the standard ground time, for the calculation of MOG 
consumption. Since the modified repair time 1s being used in the calculation of the MOG 
consumption, the MOG efficiency factor needs to be set to one. Now the model takes 
into account aircraft reliability and the analysts can examine the effect of aircraft 


reliability on the mix of aircraft used and the re-routing of aircraft. 
3. Stochastic Optimization Model 


The next step is to solve the stochastic optimization model with a larger data set. 
This thesis examined a 30 day scenario, but during this time span the airlift infrastructure 
was MOG constrained only for a period of eight days. The effects of the stochastic 


optimization should be more evident if a larger problem is used for analysis. 
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APPENDIX A. THROUGHPUT II MODEL 


The following is a brief summary of Throughput II (Morton ez. a/., 1995). 


A. INDICES 
u indexes units, e.g., 82nd Airborne 
a indexes aircraft types, e.g., C5, C141 
tt’ index time periods 
b indexes all airfields (origins, enroutes and destinations) 
i indexes origin airfields 
k indexes destination airfields 
r indexes routes 


B. INDEX SETS 


1. Airfield Index Sets 


B set of available airfields 
IcB origin airfields 
KcB destination airfields 


Zs Aircraft Index Sets 


A set of available aircraft types 

A buieCA aircraft capable of hauling bulk size cargo 
AoverGA py aircraft capable of hauling over-sized cargo 
AoutCAover aircraft capable of hauling out-sized cargo 


Bulk cargo is palletized on 88 x 108 inch platforms, which can fit on any plane. 


Over-sized cargo is non-palletized rolling stock; it is larger than bulk cargo and can fit on 
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a C141, C5 or C17. Out-sized cargo is very large non-palletized cargo that can fit into a 


C5 or C17 but not a C141. 
3: Route Index Sets 


R set of available routes 


R,cR permissible routes for aircraft a 

RapCRa permissible routes for aircraft a that use airfield 5 

Rock, permissible routes for aircraft a that have origin i and destination k 
DRicR delivery routes that originate from origin i 

RR,CR recovery routes that originate from destination k 


4. Time Index Set 


T set of time periods 

de aa possible launch times of sorties for unit w using aircraft a and route r 

The set 7,,,, covers the allowed time window for unit z, which starts on the unit’s 
available-to-load date and ends on the unit’s required delivery date, plus some extra time 


up to the maximum allowed lateness for the unit. 


C. GIVEN DATA 


A Movement Requirements Data 


MovePAX, it Troop movement requirement for unit uw from origin 7 to 
destination k 


Move UE viz Equipment movement requirement in short tons (stons) for unit 
from origin / to destination k | 


ProBulk, Proportion of unit w cargo that is bulk-sized 


ProOver,, Proportion of unit uv cargo that is over-sized 
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ProOut, Proportion of unit w cargo that is out-sized 
2. Penalty Data 


LatePenUE, — Lateness penalty (per ston per day) for unit u equipment 
LatePenPAX, Lateness penalty (per soldier per day) for unit u troops 
NoGoPenUE,, Non-delivery penalty (per ston) for unit uw equipment 
NoGoPenPAX, Non-delivery penalty (per soldier) for unit uv troops 
MaxLate Maximum allowed lateness (in days) for delivery 


Preserve a Penalty (small artificial cost) for keeping aircraft a in mobility 
system at time / 


3 Cargo Data 


UESQFt, Average cargo floor space (in sq. ft.) per ston of unit u 
equipment 

PAXW,, Average weight of a unit w soldier inclusive of personal 
equipment 


4, Aircraft Data 


Supplyar Number of aircraft of type a that become available at time ¢ 

MaxPAX, Maximum troop carriage capacity of aircraft a 

PAXSQF lug Average cargo space (in sq. ft.) consumed by a unit wu soldier fro 
aircraft a 

ACSI ta Cargo floor space (in sq. ft.) of aircraft a 

Loadl ff, Gare space loading efficiency (<1) for aircraft a. This accounts 


for the fact that it is not possible in practice to fully utilize the 
cargo space. 
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URate, 





Established utilization rate (flying hours per aircraft per day) for 
aircraft a 


5. Airfield Data 


MOGCap, 


MOGReq op 


MOGEff 


Aircraft capacity (in narrow-body equivalents) at airfield 5 


Conversion factor to narrow-body equivalents for one aircraft of 
type a at airfield 4 


MOG efficiency factor (<1), to account for the fact that it 1s 
impossible to fully utilize available MOG capacity due to 
randomness of ground times 


6. Aircraft Route Performance Data 


MaxLoad,, 


GTime aby 


DTime apy 


FltTime gy 


CTime a, 


DaysLate yar 


Maximum payload (in stons) for aircraft a flying route r 


Aircraft ground time (due to onload or offload of cargo, 
refueling, maintenance, etc.) needed for aircraft a at airfield 5 on 
route r 


Cumulative time (flight time plus ground time) taken by aircraft 
a to reach airfield 6 along route r 


Total flying hours consumed by aircraft a on route r 


Cumulative time (flight time plus ground time) taken by aircraft 
a on route r 


Number of days late unit w’s requirement would be if delivered 
by aircraft a via route r with mission start time ¢ 


DECISION VARIABLES 


1. Sortie Variables 


7 
Anar 


Number of aircraft a that airlift unit w via route r with mission 
start time during period / 
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om Number of aircraft a that recover from a destination airfield via 
route r with start time during period ¢ 


2. Aircraft Allocation and De-allocation Variables 


Allot sit Number of aircraft a that become available at time ¢ that are 
allocated to origin i 


Release gi, Number of aircraft a available at origin / in time ¢ that are not 
scheduled for any flights from then on 


3. Aircraft Inventory Variables 


Haj Number of aircraft a inventoried at origin / at time f 
AP x1 Number of aircraft a inventoried at destination & at time / 
NPlanes x Number of aircraft a in the air mobility system at time / 


4. Airlift Quantity Variables 


Tons UE yart Total stons of unit equipment airlifted by aircraft a via route r 
with mission start time during period ¢ 


TPAX ia Total number of unit u troops airlifted by aircraft a via route r 
with mission start time during period f 


=r Elastic (Nondelivery) Variables 


UENoGo xix Total stons of unit « equipment with origin / and destination k 
that is not delivered in the prescribed time frame 


PAXNoGo, i, Number of unit wu troops with origin / and destination & who are 
not delivered in the prescribed time frame 
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E. OBJECTIVE 


minimize 


»y Dy »y »y LatePenUE ,* DaysLate,.. * TonsUE ae (A.1) 


ua reR, tel, 


+) YY Y LatePenPAX,, *DaysLate jay, * TPAX sp 


ua réeR,teT,, 


+) > > (NoGoPen UE,,* UENoGo,, + NoGoPenPAX ,* PAXNoG0,, ) 
ui ok 


+) > Preserve ,,* NPlanes,, 
a ft 


The objective function minimizes the total weighted penalties incurred for late 
deliveries and non-deliveries. The model’s secondary objective is to choose a feasible 


| solution that maximizes unused aircraft. 


F. CONSTRAINTS 


> > » fons oj ane zs UENOG0 yy = Move Uk, uik ? (A.2) 
acd,y TER gy (Ele 
V u,i,k: MoveUE,,, > 9 
»y »y Sy TonsUE,,, + UENoGo,, 2 ProOut,* MoveUE,,,., (A.3) 
acd, TER tele 
Vu,i,k: MoveVE,, >9 
(A.4) 


» »y SY) TOnsVE yan + UE-NoGo,, 2 


aes FER (ET ey 


(ProOver, + ProOut,}* MoveUE,,, u,i,k: MoveVE,, >0 
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>: > »y TPAX __+PaxNoGo_, = MovePAX , , 


a reR iy tel 


V u,i,k. MovePAX,,, >0 


uik 


> es + He, , + Release, 


r reDR; 


H 


ai.t—] 


+ Alloty+Y Yay, Wait 


reR,, ¢°+[CTime,, |= 


»y Ls + HP, = cles i+) > 2X uart' ? Va,k,t 


reRR, u rer, 
al CTae, j=? 


1») Allot,, S a Y Sul, V al 


NPlanes., -y 2 Allot ji — “ 2 Release ji’, Va,t 


¥ YY Kew * Xu + DD Kant Ye + DY Hay 
i t's] 


reR, U=l u rek t= 
f 
+>) AP 32S y NPlanes,,; Val 
kos =] 
where 


. ae, if t<t<r’+CTime, -1 


wm |Clime,, if 120+ CTime,, -1 
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(A.5) 


(A.6) 


(A.7) 


(A.8) 


(A.9) 


(A.10) 








TPAX yan <= MaxPAX ,* X jan» Vu,a rere, (A.11) 


TonsUE. ., + PAXWt* TPAX ,,, < MaxLoad ,* X jan; (A.12) 
Vu,ayr,titeT 


“ar 


PAXSgFt,,* TPAX ,_, + UESqFt,* TonsUE,,, < | (A.13) 
ACSgFt.,* LoadE ff ,* X jan VuartiteT, 
YY YS FuTime,,* Xn + YY FltTime,, *¥ 44, S (A.14) 
u rek tel, reR, ¢ 
¥' URate, * NPlanes ,, Va 
YY VY (MOGRe|qy* Gime, | 24)* X sa (A.15) 


ua reéeR, vel, 
°  t'4+[DTime,, |=¢ 


+) » »y ( MOGReq,, *GTime,,, | 24)*Y,,. 


a réR, t’+[DTine,,, |=t 


< MOGLff * MOGCap, , V b,t 


A.2 Demand satisfaction constraints for all classes of cargo 
A.3 Demand satisfaction constraints for out-sized cargo 
A.4 Demand satisfaction constraint for over-sized cargo 
A.5 Demand satisfaction for troops 

A.6 Aircraft balance constraints at origin airfields 


A.7 Aircraft balance constraints at destination airfields 
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A.8 Aircraft balance constraints for allocations to origins 

A.9 Aircraft balance constraints accounting for allocations and releases 
A.10 Cumulative aircraft balance constraints 

A.11 Troop carriage capacity constraints 

A.12 Maximum payload constraints 

A.13 Cargo floor space constraints 

A.14 Aircraft utilization constraints 


A.15 Airfield MOG constraints 
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APPENDIX B. AIRCRAFT BREAK AND FIX RATES FOR 1994 


| BREAKS LANDINGS 








“7085 
169 
977 
1138 
~n7__| tad 


155 97) 15.96 
190 1045 18.18 
Ste 
1167. | ‘13.97 
1014 14.79 
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ea bag 

















ie) | | | | OO 
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FIX FIX FIX FIX 
ie RATE PIX | RATE eve | RATE aa | RATE 
0-4 4-8 8-12 | | 12-16 
| | 04 ] | 4-8 8-12 12-16 
| HRS | HRS HRS | HRS HRS HRS | HRS | HRS 
C-5 | | | | L 
8 | 195 | 8 | 19. 6 146 | 4 9.8 
FEB | 8 | 17.4 6 13.0 6 13.0 6 13.0 
MAR | 14 | 26.4 17 | 321 | 8 
APR | 16 | 25.0 18 21) 5 | 78 | 11 17.2 
MAY 15 | 29.4 11 1 21.6 3 59 5 9.8 
JUN | 11 | 13.3 16 | 19.3 11. | 133 9 10.8 
JUL 15. | 25.0 15 25.0 6 100 | 7 
AUG | 14 | 264 9 17.0 7 | 13.2 6 11.3 
SEP 10 | 27.0 10 | 27.0 4 10.8 | 4 10.8 
OCT | 17 | 304 | 11 19.6 7 12.5 5 8.9 
NOV | 15 | 224 11 | 164) 8 11.9 | 11 =| 164 
DEC | 5 oe 139 | 9 250 | 8 222 | 5 | 139 
C-141 an | i 
JAN 50 | 39.7 21.4 a 143 | 11 8.7 
FEB | 40 | 31.5 21.3 29.1 7 55 
MAR 77, —-|~ «41.0 re | 22.3 [ 15.4 | 15 8.0 
APR | 65 | 425 | 49 32.0 18 11.8 13 8.5 
MAY | 65 (| 455 43 | 30.1 15 105 | 10 | 7.0 
SUN | 74 | 46.0 4427.3 20. | 124 9 5.6 
JUL | 60 | 420 | 27 18.9 16 7 191.2 | 15 | 10.5 
Ba eh ree eee eer Cn Ae Le JOO re ene eR 
_SEP | 65 ae St 20.0 7 19 | 1230) 14 90 
OCT | 63 | 332 52 27.4 33. | 174 | 15 | 79 
NOV | 61 | 37.4 a 15.3 17, | 104 | 15 9.2 
DEC | 50 |: 333 | 32 21.3 22 147 20 | 13.3 
C-17 | ee eee eee 4 
| 
Se OO 
eee i Senerren ene heer Sere eee 
MAR | 3 —|_s37.5 oO | 00 OF Oe 25 
_APR | 2 | 667 | 0 ee a OE 
MAY | 2) 286 | 2 8G CO CO 
_ SUN) 1 33.3 ot 33.3 | = ee ee O 
IUD OO hn 0 sf OOO 0 08 
AUG O00 25.0 | Ol | CO 
OEP 0) 00 2) 00k OOO 
BOG ole a Vin OO Dh th OO. ie Oe aan 0. 0, 10 
NOV 4 40.007 oT 100 | 7 10 | 0 
DEC 3 37.5 2 2.0 | oO | oO | 0 0 
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APPENDIX C. BENDERS DECOMPOSITION CALCULATIONS 


This appendix calculates the cut gradient and intercept expressions used in 
Benders algorithm for solving the stochastic optimization model. As shown in Chapter 


IV, Section C, the algorithm will add the following cuts at each iteration 
-G,x+0, 22, Vt (C.1) 


Since x represents the first-stage decision variables used in the Throughput II model, 


constraint (C.1) is equivalent to - 
SVS Y Gt Xen + VY Gh * Lo +6, > £5, V bt. (C.2) 
ueaopr ft aor #£ 


To calculate expressions for G/,, and g,,, the first step is to take the dual of the recourse 


function. In this case, the primal is 


h,,(x,@) = min MOGPEN,,* Rf (C.3) 


br 


subject to 
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1) SY YD [Moe Reg,,*(GTime,, + MRTimeS, )/24]* X sar 


“easier ft 


abt 


+ LEE [MOGReg+( Time, + MiimeS,)/24] Yu 
ar ¢ 


— Ry S$ MOGCap, 


(2) R, 20. 
The dual of the recourse function is 
h,,(x,@) = max (C.4) 


~ MOGCAP, 
As +> YY | MOG Reqay*(GTime,y, + MRTimeS, )/ 24}* X joy 


abt 


abt 


tL YY [MOG Req,*(GLime yy, + MRTime’,) 24] Yay 
aor ¢ 
subject to 


(1) a® < MOGPEN,, 


(2) 2° >0. 
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From the L-Shaped Algorithm for stochastic programming (Morton,1994), we know 
Gi=E,n°B® and gi=E,n°d® where 1° denotes the optimal dual vector. 


Application to this model yields the following expression the cut gradient 
Gj, = E,,15,B° = E,,A3, MOG Req,,*(GTime,, + MRTime%, )/ 24 (C.5) 


and the following for the cut intercept 
24, =Eu,d° = E,,A%,(—-MOGCap, ). (C.6) 


The solution to the dual problem will depend on two possible conditions that can exist. 


In the first case, if there is no MOG violation at airfield 5, time period ¢ and scenario @, 
the corresponding recourse variable R;? will equal zero. In this case, the corresponding 
dual variable A}, is equal to zero. On the other hand, if there is a MOG violation at 
airfield 5, time period / and scenario @, the corresponding recourse variable R,, will be 
greater than zero. In this case, the corresponding dual variable A’, is equal to 
MOGPEN;,. In summary, we have 

Condition]: R, =0 A, = 0 (C.7) 


Condition Il: Ay > 0 A®. = MOGPEN,, . (C.8) 
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APPENDIX D. SIMULATION EVENT DIAGRAMS 


Appendix D contains the event diagrams for the simulation model. The following 
abbreviations are used: 

MC mission capable, aircraft not broken 

MRT modified repair time realization 


RT repair time realization. 
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Add 
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Update 
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Variables 


Notes: (1) Plane does not break, (2) Allot Orders arrive only at 
onioad airfields. 





Figure D-1: Allot Order Event Diagram. 
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Figure D-2: Release Order Event Diagram. 
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Figure D-3: X Order Event Diagram. 
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Notes: (1) Assume aircraft refueled. 


Figure D-4: Y Order Event Diagram. 
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Notes: (1) If an aircraft completes a Y order and breaks, the simulation first repairs the 
aircraft, and then ( if there exist an X order ) assigns an X order to the aircraft. 


Figure D-5: Land Event Diagram. 
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Figure D-6: Enter Service Diagram. 
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Figure D-7: Ready-1F and Ready-1P Event Diagrams. 
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Figure D-8: Ready-2F and Ready-2P Event Diagrams. 
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Figure D-9: Fill Empty Slot Diagram. 
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Figure D-10: Ready-Hold Event Diagram. 
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Figure D-11: Transition Event Diagram. 
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APPENDIX E. MOG EFFICIENCY FACTOR - LANDINGS 
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APPENDIX F. MOG EFFICIENCY FACTOR - STONS 
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